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summary 

Detailed analysis has been made of the ingestion and dispersion of 
engine exhaust products into the trailing vortex system of supersonic 
aircraft flying in the stratosphere. The- rate of mixing between the 
supersonic- jet and the co- flowing supersonic stream was found- to be an 
order of magnitude less than would be expected on -the basis of subsonic 
eddy-viscosity results. The length of the potential core was 66 nozzle 
exit radii so that the exhaust gases remain at elevated temperatures and 
concentrations over much longer distances than previously estimated. 
Ingestion started at the end of the potential core and all hot gas from 
the engine was ingested into the trailing vortex within two core lengths. 
The vortex inhibits entrainment and thermally insulates the hot gases in 
the core from the stratosphere. The final buoyancy phase started only 
ten seconds after leaving the trailing edge. The temperature rise in the 
wake of the supersonic aircraft is much greater than that for a subsonic 
transport since temperature rise was found to vary directly as speed 
squared and inversely as aspect ratio. T'»e buoyancy theory incorporates 
an entrainment coefficient which is a function of local turbulence condi- 
tions. Comparison between the buoyancy calculations for the supersonic 
case with nondimen sionalized subsonic aircraft contrail data on wake 
spreading showed good agreement. Velocity and temperature profiles have 
been specified at various!- stages of the wake, and the analysis in this 
report can be used to predict variations of concentrations of species such 
as NO under conditions of chemical reaction. 
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A STUDY OF INGESTION AND DISPERSION OF ENGINE 
EXHAUST PRODUCTS IN TRAILING VORTEX SYSTEMS 

By Jack N. Nielsen, Stephen S. Stahara, 
and James P. Woolley 
Nielsen Engineering & Research, Inc. 

1.0 INTRODUCTION 

This document represents the -final report on an analytical study of 
the ingestion and dispersion of engine exhaust products in trailing vortex 
systems under Contract NAS2-3337. The emphasis of the work is on super- 
sonic flight in the stratosphere although the analytical methods developed 
herein are for the most part not limited to those conditions. 

The problem of what influences the exhaust gas products of supersonic 
aircraft will have on the chemical balance of the stratosphere is a subject 
of great current interest and controversy. Specifically what long term 
influences the nitrogen oxides of the exhaust will have on the ozone 
balance of the stratosphere has been discussed by Johnson and Whitten 
(ref. 1-1) and by Goldsmith, Tuck, Foot, Simmons and Newson (ref. 1-2) 
who have reached opposite conclusions. One of the questions which bears 
on the general problem is how the exhaust gases from the engines are in- 
gested into the trailing vortices of the aircraft and how they are 
ultimately dispersed. For instance, at supersonic speeds the engine 
can have long potential cores which are uniformly at the exit conditions, 
and these hot cores may influence the chemical concentration of the final 
exhaust products. Also, trailing vortices can be very long for super- 
sonic flight in the atmosphere, and if the exhaust products are entrained 
by these vortices the concentrations can remain relatively high compared 
to what they- might otherwise be. Also, a knowledge of the entrainment 
of the exhaust gas by the trailing vortices is important if flight testing 
is to be conducted to obtain concentration measurements. 

In this report a physical description will be given of the ingestion 
and dispersion of the engine air by the trailing vortices including a 
breakdown of the process into several phases. The analysis of the several 
phases is then carried out in detail. Primarily the analytical methods 
are applied to the case of a supersonic two-engine aircraft flying in the 
stratosphere. 


The report contains distinct theories for the pure jet stage, the 
hot gas ingestion stage, the equilibrium vortex stage, and the asymptotic 
buoyancy stage. These theories are contained in sections 4 to 8. For the 
reader who is not concerned with theoretical details, it is recommended 
that he read section 3 which describes the ingestion and dispersion pro- 
cesses, and then-go to section 9 which contains the calculative example. 

2^0 — SYMBOLS 


Because this report contains a number of distinct theories, the 
symbol list is extensive.- In the following list, the symbols are separated 
in accordance with the sections in which they appear. If the same quantity 
is used in more than one section, it is listed under each section. In a 
few cases the same symbol has been used in different sections to designate 
different physical quantities where no confusion should result. 

Section 4.0 


A,B 
A ' ,B ' 
b 

* 

b 

c 

c_ 


c ! 


in 


constants in modified Crocco law 
constants in Crocco law 

W . * 

h /r n 

* * 

mixing region characteristic length, r e - r^ 


mass fraction of exhaust gases 
*/ * 

C ]/ C pj * 2 

TT 2 . 

y - l o / t y 

jet Crocco number, — M5 / U + 

2h j 


2 J 

specific heat at constant pressure 
constant having a value of either 0 or 1 

</**]*] 

* * 




E. 


E z /2im^ 


2 


downstream energy flux through any normal cross section of 
the mixing zone 

energy flux into mixing zone across upstream boundaries of 
zone 

* * 

H /Hj 

stagnation enthalpy 


jet stagnation enthalpy 

free-stream stagnation enthalpy 

constant in eddy viscosity equation 

eddy viscosity constant given by equation (4-50) 



length of potential core 

• It A 9 h 

m /m . 
a 3 

l/2v times mass flow rate per second entrained by jet up 
to a given axial position 

1/2tt times mass flow rate at jet exit 

Mach number on jet centerline 

Mach number on streamtube containing jet exit flow 
Mach number of jet in potential core region 
free-stream Mach number 

average mean squared velocity taking into account turbulent 
fluctuations 




radial distance from jet centerline 

radius of streamtube containing jet exit mass flow 

radius of entire jet 

radius of potential core 

exit radius of nozzle 



static temperature 

mass average static temperature over jet cross section 
static temperature on nozzle centerline 
temperature of streamtube comtaining jet exit flow 
static temperature in potential core 


ii 

I! 

II 


stagnation temperature 

stagnation temperature on jet centerline 


stagnation temperature on streamtube containing jet exit 
flow 

stagnation temperature of potential core 


U . 
min 


free-stream stagnation temperature 


static temperature in free stream 

time average mean square axial turbulent velocity fluctuation 

★ * 

u /u. 

(U - U ro )/(U C - U o ) , normalized velocity 


Ce /V j 


w w 

V"j 

* s it 

V°j 

U*/U* - 1 

* , * 
W u j 

U*. /U* 
min j 

* „ * 


axial velocity component 

mass average axial velocity over jet cross section 
axial velocity on jet centerline 

axial velocity on streamtube containing jet exit flow 
axial velocity of potential core 
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maximum axial velocity at given jet cross section 
minimum axial velocity at given jet cross section 
free-stream axial velocity 

time average mean square radial turbulent velocity 
fluctuation 

★ 

radial velocity divided by 

* ^ * 

W /Wj 

molecular weight 

molecular weight of exhaust gases 


molecular weight of free stream 


* * 

2 /r n 


axial distance behind nozzle exit 
increment in z 


P<-T- * 


CO o 


ratio of specific heats 

* * ft 


VV> 


turbulent eddy viscosity ■ e 


u 


value of e* on streamtube containing jet exit flow 


eddy thermal conductivity coefficient 

eddy diffusion coefficient 
value of e* on streamtube "s" 
eddy momentum diffusivity (viscosity) ■ e 
(r - r i )/(r Q - r^ 


(1 + Uj/(1 - U J 


density 

density on jet centerline 


density on streamtube containing initial jet exit flow 
density in potential region 
density on streamtube "s" 
free-stream density 

turbulent shearing stress associated with radial gradient 
of axial velocity 

value of t* on streamtube initially containing jet exit 
flow 

Sections 5.0 and 6.0 ; Additional Symbols 


D 


-O 


BY 


Dx 


? 

Dy 


wing span 

crossflow drag coefficient per unit length 

* * * 

unit vectors in (x ,y ,z ) direction 

nondimens ional net buoyant force per unit length, 

normalized by m^U*/r* 

] 3 p 

* 

nondimens ional drag force per unit length in x direction 
normalized by mjU^/r n 

nondimens ional drjg^fojce per unit length in y* direction 
normalized by ihjUj/r n 

* 

nondimens ional vortex force gej unit length in x 
direction, normalized by m\vJ../r* 

* 

nondimensional vortex force.gej ugit length m y 
direction, normalized hy m^U^/r n 

acceleration due to gravity 
total trailing vortex strength, 


'-’*’***• 




S ri vfir-rr 


K 


k 


in 




. * 

m. 

N 

P 

Pc 


th 

nondimensional circulation of k vortex * 1/2?', 

normalized by U*r* K 

J n 

nondimensional mass flow rate of ambient air ingested by 
jet up to a given axial location, normalized by itk 

nondimensional total jet mags flow rate up to a given axial 
location, normalized by 

initial jet mass flow rate 
number of external vortex pairs 

* 

nondimensional pressure, normalized by * 

★ 

nondimensional free-stream pressure, normalized by ij^j 


r vk 

r vmk 


‘ vmo 


“ VO 


T j 


U, 


U 


nondimensional radius of jet outer boundary normalized 
by r* 

nondimensional radial distance from jet centerline to k 

vortex, normalized by r 
* n 


nondimensional radial distance from jet centerline^to 
symmetric image of k th vortex, normalized by r n 

nondimensional radial distance from jet centerline to left 
wing panel vortex centroid, normalized by r n 

nondimensional radial distance from jet centejline to right 
wing panel vortex centroid, normalized by r n 

nondimensional time, normalized by r n / u j 

nondimensional mass-averaged jet temperature, normalized 
by T* 

★ 

nondimensional free-stream temperature, normalized by 1\ 

nondimensional mass-averaged jet velocity, normalized by U 

* 


nondimensional free-stream velocity, normalized by 

nondimensional lateral velocity in x direction relative 
to moving jet centerline, normalized by tb 


V r 

nondimensional radial flow velocity at outer boundary of 
jet, normalized by Uj 

V vx 

nondimensional velocity in ^x direction induced by external 
vortices, normalized by Uj 

V vy 

nondimensional velocity in *y direction induced by external 
vortices, normalized by Uj 

7 i 

nondimensional combination of free— stream and vortex— induced 
velocities normal to jet centerline, normalized by Uj 

w 

nondimensional velocity in y direction relative to moving 
jet centerline, normalized by Uj 

w 

nondimensional complex velocity potential, normalized by 

U*r* 

J n 

nondimensional aircraft-fixed inertial Cartesian coordinate 
system with z axis directed rearward, x axis directed 
to the right facing forward, and y axis directed verti- 
cally upward; normalized by r* 

(x,y,z) 


nondimensional crossflow coordinate system fixed to moving 
jet centerline; normalized by r* 

(x c» y c ) 

nondimensional location of jet^centerline in inertial (x,y) 
coordinates, normalized by r n 

(x vk’ y vk ) 

nondimensional location of k th right wing-panel vortex in 
inertial (x,y) coordinates, normalized by r* 

^ x vmk ,y vmk^ 

nondimensional location of left wing-panel image of k^h 
vortex in inertial (x,y) coordinates, normalized by r R 

^ x vmo ,y vmo^ 

nondimensional location of left wing-panel vortex centroid 
in inertial (x,y) coordinates, normalized by r n 

^vo'W 

nondimensional location of right wing-panel vortex centroid 
in inertial (x,y) coordinates, normalized by r n 

( x vo » y vO^ 

location of vortex centroid in dimensional inertial 
coordinates 

Z i 

nondimensional complex variable in crossflow plane, 

z i " x i + iy i 

z vk 

nondimensional complex coordinates in crossflow plane of 
right wing panel image of k^ vortex, Z vk - * vk + iy vk 
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nondimens ional complex coordinates in crossflow plane of 
left wing panel image of kth vortex, - * vmk + vm k 

nondimensional complex coordinates in crossflow plane of 
right wing-panel vortex centroid, Z vo " x vo + ^vo 

trailing vortex circulation 

polar angle in crossflow plane 

angle of jet centerline with (y,z) plane 

angle of jet centerline with (x,z) plane 

nondimensional mass-averaged jet density, normalized by p 
nondimensional free-stream density, normalized by p? 
real part of complex velocity potential W 

Section 7.0 

specific heat at constant pressure 

constants of integration 

time average stagnation enthalpy 
turbulent fluctuation in stagnation enthalpy 
laminar thermal conductivity 
turbulent thermal conductivity 

r/2v 

r c/ 27r 

value of K at position of maximum tangential velocity 
static pressure 

heat transfer rate due to turbulent fluctuations 
radial distance 

radial distance to maximum tangential velocity 
outer radius of part of vortex in solid body rotation 


r j 

r o 

R 

t 

T 


T i 

T ' 
★ 


u 

u' 

V 

v' 

v , 

w 

z 

y 


r 



p 


0 


T t 

Si 


radius at outer edge of logarithmic region 
radius where r/T 0 ■ 0.99 
gas constant for air 
time 

static temperature 
value_of T at r^ 
turbulent fluctuation in T 
free-stream static temperature 
radial velocity 

* tur bulent fluctuation in radial velocity 
tangential velocity 

turbulent fluctuation in tangential velocity 
value of v at r 

■A, i 

axial velocity 
axial distance 

ratio of specific heat at constant pressure to that at 
constant volume 

circulation, 2wrv 

total circulation of vortex 

laminar viscosity 

turbulent viscosity 

u/p 

turbulent eddy viscosity 
mass density 
laminar Prandtl number 
turbulent Prandtl number 

turbulent shearing stress, Reynolds stress 
angular velocity of vortex eye 


11 



Section 8.0 


specific heat at constant pressure 
constant of integration 
buoyant length given by equation (8-4) 
m/7TAx 

function defined by equation (8-33) 
buoyant force 
acceleration of -gravity 
function defined by equation (8-34) 
upward distance with h * 0 at t * t^ 
heat energy added to wake per unit time 
reference length given by equation (8-5) 
«Ap„c p T;(u *) 3 

airplane lift 

total mass of recirculating flow 
mass rate of entrainment 
atmospheric static pressure 

average radius of plume or wake cross section 
gas constant for air 

time 

0 

static temperature of atmosphere 

average temperature of air in wake 
axial velocity of air in plume 
free-stream velocity 

horizontal wind velocity 
vertical upward velocity of wake 


X 

Ax 

a 

e 

Y 

r 

e 

6 

Z 

H 

P» 

03 


m 

b 

c 

c 


ave 


C, 


C 


L 


m 


distance behind aircraft 

thickness of wake cross-sectional slab measured in axial 
direction 

entrainment coefficient, equations (8-1) and (8-8) 
entrainment coefficient, equation (8-1) 
ratio of specific heats 
vortex circulation 

rate of entrainment of ambient air 
inclination of plume axis to horizontal 
gVR g T* 

mass average density of recirculating flow 
atmospheric mass density 

Brunt-Vaisalla frequency, equation (8-3) 

Subscript 

initial value of a quantity 
Section 9.0 

wing aspect ratio 
wing span 

mass fraction of jet exhaust gases 

mass weighted average of c across jet cross section 

specific heat at constant pressure for air 

airplane drag coefficient 

airplane lift coefficient 

see equation (4-65) 
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see equation (4-65) 

increase in height during buoyant phase 
vr 

circulation parameter, r Q /27r 


airplane lift 

ratio of mass of entrained air to mass of exhaust gases 


Mach number on jet centerline 
Mach number in jet exit 
free-stream Mach number 

radial distance from center of trailing vortex 

*. * 

V r n 

* * 

r/r 
e n 

* * 

V r n 

radius of center to the 99-percent circulation position 


jet outer radius 

radius of streamtube originating at jet exit 
radius of potential core 
radius of jet exit 

radius to point of vortex maximum tangential velocity 

equivalent radius of recirculation region during buoyant 
phase 

initial equivalent radius of recirculation region at 
beginning of buoyant phase 

time after airplane passage 

characteristic time, equation (9-2) 




tangential velocity in trailing vortex 
vortex maximum tangential velocity 

* \ 

airplane speed (numerically equal to U w ) 
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x,y 


x v ,Y v 


i 

2* 

a. 

r o 

’i 

f^OO 


lateral and vertical coordinates, respectively, with the 
origin at the trailing-edge location on the vertical plane 
of symmetry of the airplane 

lateral and vertical coordinates with origin at centerline 
of trailing vortex 

z*/r* 

' n 

value of z for start of jet break-up 

axial distance behind wing trailing edge 

entrainment coefficient 

trailing vortex strength 

efficiency of conversion of engine heat into propulsion 
work 

free-stream air mass density 


initial value of at beginning of buoyant phase 

average mass density of fluid in recirculation region 
initial value of "p at beginning of buoyant phase 


3.0 PHYSICAL DESCRIPTION AND MODELS OF INGESTION 
AND DISPESEION PROCESSES 

3.1 Preliminary Considerations 

The exhaust products of an engine contain many components including 
water, nitrogen oxides, and particulate matter. Wake concentrations at 
the exhaust of high altitude aircraft are expected (ref. 3-1) to be 
90,000 ppm HgO vapor, 500-5000 ppm CO, 50-200 ppm NO and 1-10 ppm S0 a . 
Ambient levels at 20 kilometers are of the order of 3 ppm H a O vapor, 
0.01-0.20 ppm CO, 0.001-0.003 ppm NO, 0.0001-0.001 ppm S0 a and 2 ppm 0 9 . 

In the analysis of the problem we will neglect the particulate matter. It 
is generally sufficiently small that it will follow the streamlines, and in 
any event presents such a minor mass fraction that it should not influence 
the gas flow significantly. 
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Another assumption that we will generally use is that the temperature 
and concentration profiles in the jets are similar. Low speed theory and 
experiment show that the transfer of heat and mass are similar but 
both differ from the transfer of momentum. The radial spread of tempera- 
ture and concentration is slightly faster than that of velocity. 

Wilder and Hindersinn (ref. 3-2) have shown by measurement that the 
spread of a C0 2 jet exhausting into a supersonic airstream was in good 
agreement with jet spread data, determined on the basis of total tempera- 
ture measurements. They also showed- that gas dynamic expansions and- 
compressions within the jet caused no essential change in the nature of 
the jet spread due to turbulent mixing. The rate of spread due to mixing 
was found to be small. Data on free jets generally show similar concen- 
tration and temperature profiles (ref. 3-3, pp. 98-102). If the transfer 
of heat and mass is dominated by turbulent diffusion and the various gas 
components of the engine exhaust remain in the same relative proportions, 
it is reasonable to take all concentration profiles similar to the 
temperature profiles as a first approximation. Accordingly, we are faced 
with determining the temperature field of the flow. 

Consider the overall flow process to consist of a number of stages. 
The first stage is the jet stage in which the jet issues from the engine 
and is relatively unaffected by the presence of the wing vortex wake for 
a distance downstream of the aircraft. In the next stage the trailing 
vortices have sufficient influence to break up the jets and entrain some 
or all of the engine exhaust. We will designate this stage as the hot 
gas ingestion stage. After the ingestion stage there exists the 
equilibrium vortex stage in which the vortex pair act together and move 
with a regime of recirculating flow which is closely elliptical in cross 
section. In the final asymptotic stage the buoyancy force acting on the 
recirculating flow together with turbulent entrainment of atmospheric 
air combine to produce the asymptotic stage of the motion. 

3.2 Jet Stage 

The initial conditions for the study are taken at the point where 
the engine exhaust exits to the atmosphere. For the example calculations 
it is assumed that the jet pressure is equal to the atmospheric pressure 
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at the exit. This Is achieved by means of engine flaps. If pressure 
equalization takes j-lace external to the nozzle, a special calculation 
may be required to determine the initial conditions corresponding to 
atmospheric pressure. In addition to pressure equalization at the exit, 
we have also assumed that the velocity profile is uniform. 

At the exit the flow is described in terms of a static pressure, 

Mach number and stagnation temperature. The jet exhausts into a 
mainstream at atmospheric pressure and temperature and airplane Mach 
number. The jet Mach number may be less than that of the free stream 
because of the high speed of sound in the jet. 

Consider first the case where the jet axis is aligned in the flight 
direction. The supersonic jet in a coflowing supersonic stream now 
entrains free-stream air and slows down. For a certain distance behind 
the aircraft, the potential core exists and beyond this distance the "far 
field" occurs. The velocity and temperature profiles are governed by 
turbulent mixing between jet and free stream and on the stagnation 
temperature ratio of the jet and free stream. No reliable method of 
calculating the properties of the jet and based on data for hot super- 
sonic jet in coflowing supersonic streams was found in the literature. 
Accordingly, one has been developed, as described in section 4. It turns 
out that the entrainment rates are much less than expected. 

It is only at certain angles of attack that the jet thrust axis is 
aligned with the free-stream direction. For other angles of attack 
there is thus the additional complication of a jet obliquely inclined to 
the free-stream flow. The analysis of Wooler (ref. 3-4) for oblique 
subsonic jets has been adapted to treat this effect by means of a 
crossflow drag coefficient. 

There is also an effect of the trailing vortices on the jet as it 
leaves the aircraft. One effect is that the crossflow velocity normal 
to the jet is reduced because of the induced downwash between the vortices 
so that the crossflow drag on the jet is reduced. It is also assumed 
that the jet before it disintegrates acts as a solid boundary for the 
purpose of imaging the trailing vortices. There is thus an additional- 
component of force acting on the jet. Furthermore the motion of the 
vortices induced by the presence of the jets cause an additional force 
component acting on the jet. Finally there is the buoyant force acting 
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on the jet by virtue of its higher temperature than the surroundings. 

As a result of all these forces, the jet tends to straighten out parallel 
to the trailing vortices after a certain distance of travel behind the 
wing trailing edge. At some point near this position the hot gas 
ingestion stage commences, and the jet is allowed to disintegrate and 
become entrained by the trailing vortices. 

3.3 Hot-Gas Ingestion Stage 

No data have been found to illustrate how the jet from an engine 
breaks up and is entrained into a trailing vortex. In order to obtain 
some insight into the question, a test was made in the 40- by 80-Foot 
Wind Tunnel using the 6- foot span Convair 990 wing model designed and 
built for a trailing- vortex investigation. A trailing vortex was made 
visible by placing white smoke generators near the wing tip. At about 
the 69-percent semispan location slightly under the wing, a yellow 
smoke generator was used to simulate the engine jet. Smoke was forced 
out of the engine location at a low tunnel speed (30 ft/sec, approximately) 
so that observations could be made by investigators in the wind tunnel. 

The yellow smoke moved downward under the vortex, and reached a point 
directly under the vortex at about 30 feet downstream of the wing. As 
it moved past this point, it rapidly broke up and became ingested by the 
trailing vortex. The engine smoke thus maintained its identity for 
about five wing spans. 

In calculating the jet break-up and its ingestion into the trailing 
vortex we have considered the jet no longer to be under the influence 
of crossflow drag nor to contain an image vortex system. In fact, the 
various regions of the jet are now considered to be independent and to 
be acted on only by the direct effect of the trailing vortex pair and 
by buoyancy. The pieces of the jet considered as mass points then tend 
either to spiral into the vortex or escape from the recirculation region. 

During this process the initial temperatures of the various 
pieces are obtained from the calculated temperature profiles in the jet. 

As they are ingested, it is assumed that their temperature time history 
is the same as if the jet had continued to develop. The turbulent mixing 
processes in the two cases are probably different, so that the assumption 
is only approximate. However, if the hot gas is rapidly ingested, then 
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this assumption should be sufficiently accurate since the gas will not 
cool much during ingestion. 

The above calculation should show how much of the exhaust gas is 
ingested. The amount of exhaust gas can be much greater than the amount 
of air in the core of the undisturbed vortex. Accordingly, hot gas 
ingestion will change the structure of the trailing vortex. All the hot 
air tends to be concentrated in the center of the vortex. Some effect 
of hot air_ingestion on vortex stability might be expected. 

3.4 Equilibrium Vortex Stage 

In the equilibrium vortex stage the ingested exhaust gases become 
part of the recirculating flow around the vortices. The general flow 
picture is shown in the following sketch as taken from reference 3-5. 



The recirculating flow near the vortices will be very similar to that 
for a single trailing vortex since the interference of one vortex on 
the other is small under these conditions. The whole recirculation 
region mass moves upward or downward in the atmosphere as a whole. The 
calculation of reference 3-5 of the rate at which vorticity diffuses 
from the recirculation, region into a surrounding quiescient atmosphere 
or is annihilated along the vertical axis of symmetry shows a decrease 


20 


in the vortex strength in the ratio of 1/e in about 80 minutes. While 
this case is for a subsonic transport, it should also give a good order 
of magnitude estimate for a supersonic aircraft. 

In the case of mass diffusion of the exhaust gases from the 
recirculation region, mass can be detrained from the recirculation region 
only across its outer boundary whereas in the case of the diffusion of 
vorticity across the vertical plane of symmetry, vorticity of opposite 
sign can annihilate each other. Considering now the analogy between 
the diffusion of vorticity and mass, we would expect the exhaust gases 
to be detrained from the recirculation region even more slowly than the 
above calculation indicates for vorticity since a means of annihilating 
vorticity was present along the vertical plane of symmetry where no such 
means is available to annihilate mass. 

The rate of escape of exhaust gases from the vortex is small during 
the equilibrium vortex stage. The vortex tends to retain vorticity and 
fluid within its confines. The vortex core is highly stable with a 
tendency to damp out turbulence resulting in very low rates of diffusion. 
The exhaust gases will concentrate in annulii surrounding the core at 
radial distances according to the relative density between the species 
and the environment. The vortex will act as a centrifuge in which gases 
will be stratified according to the radial variation in density. 

3.5 Asymptotic Buoyant Stage 

In the asymptotic stage the recirculation region is influenced 
chiefly by buoyancy. Accordingly it will tend to rise. At the same 
time its vertical motion will cause entrainment of stratospheric air. 

A detailed analysis of the asymptotic motion of the wake in a quiescient 
stratosphere is made in Section 8. 

An approximate analysis of the asymptotic stage of the wake motion 
has been made by Conti, et al. (refs. 3-6 and 3-7) considering a number 
of separate nonlinear effects. 

These authors suggested that, in addition to the buoyant rise and 
wake-induced turbulent growth, further dispersion occurs because of 
wake collapse and horizontal wind shear. According to their analysis, 
buoyant rise is the least influential of the mechanisms giving rise to 
horizontal wake growth. 
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The towing tank experiments of Fay and Overcamp (ref. 3-8) modeled 
an atmosphere without wind shear, wind generated turbulence, and strati- 
fication. The experimental results confirm the predicted growth due to 
buoyancy in the absence of these atmospheric effects. Although wake 
collapse has been observed in the case of wakes in a stratified medium, it 
has not been observed in experiments on buoyant plumes. When buoyant plumes 
level off, they become internally stratified to the same extent as the 
surrounding fluid and no further turbulent mixing occurs. It is still not 
clear whether wake collapse will play a significant role in the asymptotic 
stage, but is seems probable that buoyant forces play a more significant 
role for supersonic aircraft than for subsonic aircraft because of greater 
-:ake heating. 

The ultimate dispersion, on time scales of months or longer, of 
exhaust products will be determined by the rates of natural stratospheric 
dispersion which are governed by horizontal wind shear and stratification. 
The present study is mainly concerned with establishing the initial con- 
ditions for the dispersion by considering the movements and rates of 
dilution of the exhaust products from the engine exhaust through the 
aircraft wake. 

4.0 ANALYSIS OF HOT SUPERSONIC JET CHARACTERISTICS 

4.1 Introduction 

A method of determining the hot jet characteristics for a supersonic 
jet in a coflowing supersonic stream verified by experiment for the range 
of interest was not found in the literature. Since the details of the jet 
flow must be known to determine whether it is entrained by the trailing 
vortices, it was necessary to develop a theory for- .the jet characteristics. 
Only one set of data was found (ref. 4-1) for supersonic jets in supersonic 
coflowing streams with heating, and these data have been used to determine 
the eddy viscosity values which are the basic empirical inputs into most 
if not all jet theories. 

In the theory which follows constant pressure turbulent mixing has 
been assumed. For our example case the jet exit pressure and atmospheric 
pressure are closely matched. We have assumed that the velocity, concen- 
tration and temperature profiles at the jet exit are uniform. In the case 
of pressure equalization external to the nozzle some departure from the 
conditions of uniformity is to be expected. 
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In the ensuing analysis most of the work is carried out in terms of 
nondimens ional quantities. Such quantities as density, temperature, velo- 
city, and stagnation temperature have been nondimensionalized by the 
values of these quantities at the jet exit. All distances have been non- 
dimensionalized by the jet exit radius. Asterisks have been used to 
designate dimensional quantities and nondimensional quantities carry no 
asterisk. 


4.2 Basic Equations of the Jet 

The problem to be analysed is the turbulent mixing of an initially 
uniform axisymmetric jet exhausting into a uniform coflowing stream under 
the condition of a uniform static pressure throughout the flow field ex- 
ternal to the jet nozzle. The general features of the flow field as well 
as the nomenclature are depicted in figure 4-1. The axisymmetric flow 
field is described in terms of a radial and an axial coordinate, r* and z*, 
and the jet nozzle radius, r n » The local thermodynamic and gas-dynamic 

parameters are designated as follows: the density, p*; the stagnation 

* 

temperature, T Q ; the Mach nuntoer, M; and the velocity, U*. The free- 
stream properties are denoted by the subscript, « , while those of the 
initial uniform jet are denoted by the subscript j. The properties 
evaluated along the axis (r =0) are denoted by the subscript, c, those 
evaluated on the streamtube containing the original jet mass flow are 
denoted by the subscript, d, and a general streamtube by the subscript s. 
The outer boundary of the mixing region is designated by the radius r Q . 

The jet is considered in two parts. For some distance downstream of 
the nozzle exit a potential core exists, as shown in figure 4-1, in which 
the velocity and thermodynamic properties are uniform at the values at the 
nozzle exit. The local radius of the potential core is designated r^. 

The axial distance over which the potential core exists will also be termed 
the near field. Downstream of the potential core the velocity on the jet 
centerline U c progressively decreases. This region of the jet is '..erm'-d 
the far field. 

The integral form of the governing equations for conservation of mass 
and momentum along a streamtube, s, are: 
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Mass 


dz 


r 

/ s * * * * 

p U r dr =0 


(4-1) 


Momentum 


d r s *- * _ * * 

-A. p (U ) 2 r dr 

dz J 


***** 
r g (z ) t (z ,r g ) 


(4-2) 


where 


r (z ) is the radius of a given streamtube, s, at any axial 


position, z . 

T*(z*,r*) is the shear stress on the streamtube, s, at the axial 
s * 
position, z . 

* i 


* * au 

p 3 e s “* 

3 3 dr 


for turbulent flow 


e is the eddy viscosity in the flow field evaluated at 
the streamtube, s. 

Integration of equation (4-1) along the streamtube beginning at the 
lip of the jet nozzle results in 


/ 


★ * 
r*(z ) 


d ' * * * * * * *n * * 

p U r dr - PjUj ~ 2 ~ - m^ 


(4-3) 


where m^ is 1/2 tt times the mass efflux from the jet nozzle. 

The total jet mass flow per unit time is proportional to the quantity 


* * 

/ r e ^ it it it * it it if it f® ★ * * * • it 

p U r dr - J P U r dr + J* p U r ur • + m a (z ) 


★ it 

?d (z > 


*/ *\ 
U ) 


(4-4) 
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of free-stream fluid entrained 


where m*(z*) is 1/2 tt times the mass flux 
in_.the_. mixing region up to the station z . 

At the outer edge of the mixing region, r 0 (z ) the shear stress, 
T*(r*(z*)), vanishes. The total jet momentum is thus proportional to the 


★ * 


quantity 


r 

/ 


* *x 

r (z ) 

-e **«*.* • * * •*.*,.,* 

p U 2 r dr = m j u j + m a ^ z ^ U «> 


(4-5) 


which indicates that the momentum inside r g at any axial station is equal 
to the original jet momentum, 


• * * * , *v o 

- irp.. (u..) 2 r n 2 

plus the momentum ingested into the mixing region from the free stream, 


* * * 
2™ a U oo- 


If we replace m using equation (4-4) 


.* r e * 
m a " J P U 


e * * * * . * 

r dr - nij 


(4-6) 


Equation (4-5) may be put into the form 


/ e * * * *, *, * 
p U (U -Uj.E.dr 


• * * 

^(Uj-Uj 


(4-7) 


for future reference. 

A similar treatment of the energy equation yields 

★ 


JL 

/ e * * * * 
p U H r 


* . * * •*„* 
dr - m^Hj + m a H w 


(4-8) 


where H* is the stagnation enthalpy of the fluid. By substitution of 
equation (4-6) into equation (4-8) one obtains 
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i 


f 


e * * * * * * •*, * 
p D (H -Hjr dr - m.. (H ^ - HJ 


(4-9) 


We note that it is unnecessary up to this point to presume that the dif- 
fusion of momentum and energy have proceeded to the same outer boundary. 
The above relations merely require that r g encompass the more extensive 

of the two, so that there is no contribution to either integral of equa- 

* * 

tions (4-7) or (4-9) for r > r Q . 

It will be convenient to normalize all quantities with respect to 
their values in the initial uniform jet since all such quantities will 
always be nonzero which may not be the case for the free stream. The 
resulting quantities will be simply those previously defined, but with 
the asterisk removed. 



(4-10a) 


for a perfect gas of uniform molecular 
weight in a constant pressure field ' “ 
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* 

U 

* 

U i 


H 


H_ 

* 

H, 


o 
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T o’ 


for a perfect gas of uniform 
specific heat 


(4-10c) 


(4-10d) 



(4-10e) 


(4-10f) 


With the assumptions of uniform molecular weight W and specific 
heat c p and using these normalizations, we can write using the ratio 
of specific heats -y * c p^ c c * 


j . 

t '■ 


fi 

j 


j 


i 

« 


i 

\ 
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‘ ipfir' 


Too - To <1 + '"T M i )/(1 + 2^ M » ) (4 “ X 

00 

o. - M^Mj (4-1 

We also normalize the entrained mass flux, 

m a - m>; (4-1 

The governing equations may now be rewritten in nondimensional form 


Mass 


?d (z) 


$ r dr - 1 


r (z) 


m a (z) 


£ r dr - 1 
T 


Momentum 


r (z) 


2 (0 

T 


Ujr dr 


1 - u 


_a r dU) s£ , dc . r a ( ‘<3 aa| 

az-J T ' dt T d (* 0 *J 3r 


For a thermally perfect gas with W = 1 the energy equation may be 


written as 


H„) t dr 


- 1 - M. 






and for a calorically perfect gas with uniform specific heat (i.e., 
c p - 1) 


r e (z) 

2 f § (T - T 0 )r dr - 1 - T c (4-19) 

J 1 O „ on 


In order to complete the formulation of the problem, it is necessary 
to introduce further information on U, T and The former two appear 

in integrated form only and reasonable approximations to their functional 
dependence on r and z should yield a reasonably good solution. The 
axial coordinate enters the problem explicitly only through equation 
(4-17) . It is thus seen that the eddy viscosity determines the axial 
scale of the solution. 


4.3 Temperature-Velocity Relationship 

Equations (4-14) , (4-15) , (4-16) , (4-17) , and (4-19) can be evaluated 
if nondimens ional similarity profiles for velocity and temperature are 
introduced into the integrals. Another approach would be to utilize a 
nondimens ional similarity profile for either temperature or velocity, and 
relate these variables through a Crocco-type relationship. If the latter 
course is chosen, then equation (4-19) is superfluous. To the extent the 
Crocco-type relationship obeys the law of the conservation of energy, 
equation (4-19) will be fulfilled. In our analysis we have chosen to use 
a Crocco-type relationship between temperature and velocity and to make a 
check of equation (4-19) with the resulting solution. Assuming unity 
Prandtl number, approximate relationships between temperature and velocity 
have been given by Crocco (ref. 4-2) and Woolley (ref. 4-3) for laminar or 
turbulent flows. 

The Crocco relation indicates that the stagnation temperature of a 
thermally and calorically perfect gas having Pr - 1 and a constant and 
uniform specific heat is a linear function of velocity. 

T* - A' + B'U* (4-20) 


The coefficients A 1 and B' may be suitably defined to satisfy the 
boundary conditions on T* and U*. 
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This relation has been well confirmed for low speed laminar flows. 

For high Mach number flows, however, it has been found to be at consi- 
derable variance with the measured temperature-velocity distributions. 

In particular, in the so-called hypersonic limit the stagnation temperature 
has been found to vary with the square rather than as a linear function of 
the velocity. 

Woolley (ref. 4-3) has shown that for turbulent flows in particular, 
the time mean static temperature for a perfect gas with constant and uni- 
form-specific heat should be a linear function- of the time mean velocity 

T* « A + BU* (4-21) 

for unity turbulent Prandtl number with the turbulent processes dominating 
the molecular ones. Both T* and U* are time average quantities and the 
coefficients A and B are evaluated for the particular boundary conditions. 
It is noted that this relation will approach the observed stagnation 
temperature— velocity relation for high Mach number flows} i.e., 

T * q 2 * .. 

T - A + BU + ■ a — (4-22) 

° 2c 

/C P 

where (q 2 ) " - u 2 * + (u' 2 )* + (v' 2 ) * -* (U 2 )* assuming (U 2 )* » (u' 2 ) 

+ * w here (u' 2 ) * and (v' 2 )* are the mean square of the longitudinal 

and transverse components of the fluctuating velocity field, respectively. 

This modified Crocco relation, equation (4—22) , was shown by Woolley 
to satisfy the time mean energy equation appropriate to turbulent flows 
with unity Prandtl number. 

Both the Crocco and the modified Crocco relation can be expressed in 
a single form valid for the entire jet for the following boundary conditions 
at the jet exit and at the frec-stream boundaries. 

Jet Exit: T G ^ - Tj » Uj " 1 

Free-Stream Boundaries: Given T 0 , T^, U^, of 

00 

which 2 are independent. 
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The stagnation temperature is of the form 


T - T 

^2 

1 - T„ 


U _U cc 
c c© 

1 - u 


DCJ^ 


1 - T, 




- DC^(U -UJ 2 

u + J i-V — u2 


(4-23) 


where 


C 2 

D 


* 

V 

jet Crocco number * 

2Hj 


- 1 


m3 


l + 


2 ”3 

T -- i 


2 M j 


0 for Crocco relation, equation (4-20) 

1 for modified Crocco relation, equation (4-22) 


U 


normalized velocity ■ (U - U ro )/(U c - U w ) 


The static temperature form is 


T - T 

0* 

1 - T 


U - U 

C oo 

1 - u 


+ (1 - D) 


1 Z 1 u2 

_2 Hi 

1 - T 


(1 - 


UJ (u c - Uj 


u 


- (1 - D) 1 2 . T j (U - U W ) 2 U 2 (4-24) 

00 

In the above expression we note that U Q ■ ■ 1.0 in the potential 

core region (i.e., for z £ £ ), but U is variable for z > £ . 

pc p 

4.4 Reduction of Equations to Similarity Form 

It has been experimentally observed that the velocity and temperature 
profiles in jets exhibit self-similar behaviors. That is, appropriate com- 
binations of the radial and axial coordinates will reduce these profiles to 
functions of a single variable to a high degree of accuracy. The resulting 
mathematical advantage afforded by such a simplification is that many of 
the integrals in the governing equations then become independent of z, and 
require calculation only once. Hence we will seek to pose the problem in 
such a form. Since we have already expressed the temperature as a function 
of velocity, we are concerned only with the velocity similarity law. 

Numerous theoretical and analytical investigations of turbulent mixing, 
Schlichting (ref. 4-4), Abramovich (ref. 4-5), Crane (ref. 4-6), Maydew 
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and Reed (ref. 4-7), have indicated that the velocity profiles for both 
incompressible and compressible flows are virtually indistinguishable. 

The functional form suggested by Abramovich for axisymmetric flows will 
be used: 

U - (1 - n 3/a ) S * 0 1 T) 1 1 (4-25) 


where 


U = normalized velocity distribution 

« (U - UJ/(U C - Uj = (U* - t£)/(U* - UJ (4-26) 

tj * similarity variable 

“ (r - r i )/(r g - r^) = (r - r^/b (4-27) 

This definition of terms in equations (4-26) and (4-27) permits the 
continuous application of equation (4-25) over the entire length of the 
jet. In the region of the potential core, U c ■ Uj ■ 1. In the far field 
beyond the potential core, U Q - U Q (z) and r i « 0 with z -Sp- 
in reducing the jet equations to similarity form we first note that 

r - b(rj + r^/b) (4-28) 

£-\ . b #- (4-29) 

dr |z ^ z 


A 
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- &| n + $s| r &|, - 


A. 

dz 


. „abV-| 

b \ dz ^ dz/di] L 


(4-30) 


and 


f f fi h 

/ g(r)dr » / g(r)dr + / g(r)dr ■ I g(r)dr 


if g(r) - r for r r^. 
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We will now rewrite equations (4-14) , (4-15) , (4-16) , and (4-19) in 
a similarity form using equations (4-25) through (4-31) . 


Mass 


and 


b* 



§ *1 dr| + 




(4- 


m a U) 


1 i 

2b s J t] dt) + 2r.jb j ^ dr] - 


(1 - r?) 


(4- 


Momentum 
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Energy 
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The conversion of equation (4—17) to similarity form is somewhat more 
complicated. The left-hand side of equation (4-17) may be expanded as 
follows: 


jL 

dz 


r Az) r r L [z) r (z) 

J f '«-■& \ / 

i lo z L ( Z) 


1 £ 

T 


r d’ 


(4 


The first term in the brackets on the right-hand side represents the 
axial momentum of the core. 


_d_ 

dz 



(*) 

ui 

T 


r dr 



. ffi 

i dz 
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since ■ 0 when U Q and T q / 1. 


The second term on the right of equation (4—35) may be written 


?d (z) 9 U? dr. 

f a U 2 . d _ d 

J T r dr “ T d r d dz 


dr. ?d (z) 


ar . r 

- *t -j? + j 


c) (U 2 /T) 


r dr (4-37) 
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The last term in the above may be written 


?d (z) 


i ± (z) 
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Noting that 


d(U 2 /T) 

dU_ 


<= 0 when r . / 0 . 


the first term on the right of equation (4-38) becomes 




r, (4-39) 
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integrating by parts and noting that 


v * Id 
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The second term on the right of equation (4-38) becomes 



(4-40) 


The right-hand side of equation (4-17) is rewritten simply as 



r n T d 



(4-41) 


^d 


Collecting equations (4-35) through (4-41) we may now rewrite equation 
(4-17) in similarity form 
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(4-42) 


where 



(4-43) 


In the application of equation (4-42) it is recalled that either 


or 


dU 

- 0 , ^ > 0 lear field) 



(far field) 
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4.5 Eddy Viscosity Determination 

The axial length sca.e through equations (4-17) or (4-42) is locally 
proportional to the eddy viscosity, e d * There is, unfortunately, no 
reliable theory and little direct information on the determination of the 
eddy viscosity for supersonic jets in supersonic co-flowing streams. Our 
approach is to correlate such data as can be found z. i the literature into 
a model for the eddy viscosity, e* and on the basis of that correlation 
estimate the eddy viscosity for the present circumstances. 

The eddy viscosity is defined by the relation 

T * s pV &L. (4-44) 

r)r* 


where 


t* = turbulent shear stress (Reynolds stress) associated with 
a point in the flow field. 

* . 

p « density at the point. 

dU*/dr* = radial gradient of the velocity at the point. 

The form of eddy viscosity used extensively for incompressible free 
shear flows is 


* * * 

Kb U mav - U . 

max mm 


(4-45) 


where 

b* * measure of the shear layer thickness 

K * nondimens ional coefficient whose actual magnitude depends 

* 

on the choice of b . 

U* ,U*. « respectively . the maximum and minimum characteristic 

max* mm r 1 * 

velocities of the shear flow. 

In the present analysis we have adopted the general functional form 
of e* given in equation (4-45) . We have taken b to be the radial 
width of the mixing region, 

b* « r* - r* (4-46) 

e l 
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and 


* * * * 

u - u . = u - u 

max mm c <« 


14-47) 


For the present analysis, any effects on the eddy viscosity rot fully 
accounted for by these two factors are contained in the parameter K. 
Thus, 


e* = Kb* (U* - U*) (4-48) 

Numerous investigators have noted that the eddy viscosity is dependent 
not only on the velocity difference across the free shear layer, but also 
on the velocity ratio, , (see, for example, Schlichting (ref. 4-4), 
Abramovich (ref. 4-5), and Sabin (4-8)). This additional dependency on 
velocity ratio appears to be expressible as a function of 

X = (1 + UJ/(1 - UJ (4-49) 

The effects of density differences in the far field in the absence of 
thermal differences was investigated by Ricou and Spalding (ref. 4-9) and 
found to be expressible as a function of the density ratio, ,j k . 

Investigations with heated jets, e.g., Corrsin and Uberoi (ref. 4-10) 
have indicated that the temperature ratio of the streams, T^ , is a cor- 
relating parameter in this case. 

Warren (ref. 4-11) found little effect of temperature on the eddy 
viscosity when the effect of Mach number was included in his study of 
single supersonic and subsonic jets in a quiescent atmosphere. It now 
appears that this effect is interchangeably expressed with the combined 
effects of the density and total temperature ratios. Similar data from 
Maydew and Reed (ref. 4-7) are in accord with this observation. 

Wilder and Hindersinn (4-12) studied the mixing of a two-dimensional 
jet in a co-flowing stream. The total temperature ratio and the Mach 
number of both streams were varied systematically in their study. The data 
are the only ones of their type found in our literature search. Their 
results are detailed in Wilder and Hindersinn (ref. 4-1) and provide 19 
data points within the following ranges: total temperature ratios, 
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0,5 T 0 0.8 ; jet Mach number, 0.9 < M . < 2.8 ; and free-stream Mach 

number, 1.25 <. <. 1.6. An eddy viscosity correlation depending solely 

on the ratio of the two Mach numbers (i.e,, independent of To and T ) 

is suggested in this study. This correlation, based on the ratio Mj/M^ , 

is clearly not applicable to the condition of a quiescent free stream where 

M_ = 0. 

00 

With this background K is assumed to be of the form 

K = K 0 f(H M ,T 0eo ,U eo ) (4-50) 


where it is assumed that 


f (Poo. T o ro » u J = P^o^ (4-51) 

K « absolute constant 
o 

An attempt was made to correlate all available data within this format 
and to evaluate K Q and the exponents a,b, and c. 

The eddy viscosity data used in the correlation study were put into 

the format indicated by the preceding equations. Because of the many 
. . . * * 

definitions of € ,b , and the analytical approaches used in the various 
investigations, it was necessary to reduce the published data to the 
present format by the procedures given in Appendix 4A. The data were all 
reduced to give the effective K defined in equation (4-50) for the given 
experimental conditions with the definition of terms prescribed in the 
preceding equations. The resulting data and sources are given in Table 4-1. 

The data of Table 4-1 were systematically used to evaluate the expo- 
nents a,b, and c of equation (4-51) . These data were primarily determined 
from near field measurements. It was found that the coefficients a,b, 
aid c were all approximately equal within the accuracy of the given data. 
The values of K may thus be correlated as a function of the single 
parameter: 

a - Hro T 0 * (4-52) 

00 

The correlated data are shown in figure 4-2 where it i3 observed that 
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(4-53) 


f 


<P»» T o » U J 
00 


a 


0.9 


appears to express the dependence of K on a very well. The free- jet 
data of Warren and the data of Wilder and Hindersinn both appear to fall 
along lines parallel to a 0,9 


The Mach number dependence noted by the preceding investigators 
appears to be reconciled by equation (4-53) . To see that this is reason- 
able we note that 


a 





X 


a 


1 + 
1 + 




X 


(4-54) 


From equation (4-54) it is seen that for Warren's free- jet case where 


M,,, = ^ “ 1* the effective K was found to decrease with increasing 

Further expansion of equation (4-54) shows that a is only a weak 
function of T 0 when expressed as a function of M and M.. It is 

00 00 j 

noted that 





M 1 V (1 + MJ) 

Combining this expression with equation (4-54) yields 
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Examination of equation (4-55) indicates that a is almost independent 
of T 0oo for small values of M^/M^ and only a function of T 0 V 2 other- 
wise. Thus’ the correlations of eddy viscosity using Mach number will show 
little effect of stagnation temperature in accord with both the results of 
Warren and of wilder and Hindersinn. 

A: can be seen from figure 4-2, the correlation of K as a function 
of a alone shows two distinct curves for M TO * 0 and M^ > 1. It was 
noted that the upper group of data (i.e., those with the higher values of 
K q ) were obtained with the external stream either quiescent or with the 
entire flow field incompressible. The lower group of data (the 19 points 
of Wilder and Hindersinn' s data) were generally taken when both of the 
streams were supersonic. Three data points from Eggers and Torrence 
(ref. 4-13) fall between the two extremes. Two of the three points are 
for > 1 but M j < 1 while for the third point, both streams were 
subsonic. It thus appears that the mixing of the jet with the co-flowing 
stream depends strongly on whether M K ■ 0 or M w > 1. It is possible 
that this result is associated with the different signs of the rate of 
change of streamtube area with axial distance for decelerating flow at 
subsonic and supersonic speeds. 

The range of a in the supersonic jet and supersonic free-stream 
data (Wilder and Hindersinn) is approximately a factor of 20, from 1.55 
to 29.5. For the case of present interest a is approximately 4. It is, 
thus, in the range of the Wilder and Hindersinn data, whereas the free-jet 
data are for a <. 1. 

The preponderence of data for conditions closest to those of present 
interest, even though these data are almost entirely from one investigation 
suggest that the eddy viscosity coefficient, K, for mixing of two super- 
sonic streams is best represented by 




K - 0.0009 a • 


(4-56) 


Therefore, this is the formula employed in the analysis. 


4.6 Computer Solution of Equations 

The equations previously developed were programmed for computer inte- 
gration in the streamwise direction to obtain the development of the flow 
field. It is noted that in the near field the knowledge of one local char- 
acteristic of the flow field, e.g., b, will permit the calculation of the 
remaining ones, e.g., r^ and r^ , through equations (4-33) and (4-31). 
Similarly, in the far field, if one knows, for example, U c , the same two 
equations permit the computation of b and r^ . Thus, it is only necessary 
to obtain the variation of one characteristic with axial position to obtain 
all. This variation is obtained by integrating equation (4-42) in an 
appropriate form. 

4 .6.1 Near field calculations .- In the near field, U c » 1 and 
dU c /dz « 0. Thus, equation (4-42) may be written as 





(4-57) 


The derivatives, dr^/db and d^/db are obtained by differentiating 
equations (4-33) and (4-31) , respectively, with respect to b . 




(U - UJ T] dq + 



(U 



(U - Ujd n 



(4-58) 


(4-59) 


The value of is determined at each station by iteration of equation 

(4-31) with the values of b and r. already computed. 



The integration of equation (4-57) was carried out by assuming a step 
size Az and determining 

db 

b(z + Az) = b(z) + Az dz (4-60) 

The initial conditions are b(0) m 0, r^(0) = r d (0) ® We then 
determine r^ by 

db ^ i 

r ± (z + Az) = r^z) + Az (4-61) 

This value of r^ was compared to that determined from equation (4-33) 
for b(z + az) . If the two values did not agree within an allowable 
tolerance (e.g., 0.1 percent), Az was halved and the computations repeated 
until the tolerance was met. The calculation of Ti d was then performed by 
iteration of equation (4-31). New derivatives, dr./db, d-r) d /db and db/dz*. 
are then computed and the cycle continued until the position r^ = 0 is 
found. This is the end of the potential core. 

4.6.2 Far field calculations .- In the far field r^ * dr^/dz = 0 and 
U c is used as the primary dependent variable in the integration of equa- 
tion (4-42) . In a procedure identical to that for the near field the 
derivatives db/dU c and dt| d /dU c are found from equations (4-33) and (4-31) . 



(4-64) 


Equation (4-64) was used to evaluate AU C at a selected Az as was 
done with b in the near field. The initial condition " * *- s 

used with the remaining values coming from the near field computations. 
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A tolerance was set on b determined from Az , equations (4-64) and (4-62) 
and from equation (4-33) as was done for r i in the near field. Upon 
satisfaction of the tolerance, other values were calculated at the new 
position. 

The new value of t| d is again computed from equation (4-31) for U 
and b at the new station. 

In the calculative procedure the energy equation is never required. 

It is, however, important that the assumed velocity-temperature relation 
satisfies this equation. This is checked at each step (Az) of the calcu- 
lation in both the near and the far fields. Equation (4-34) is put in 
the form 


or 

\ (1 - T 0 ) (1 - rf) 

Z 00 1 

+ / T ’l dT l + br i 

“L o 

Both the left and the right-hand members are calculated and printed out 
for comparison at each step of the calculation. It is seen that the left- 
hand side is related to the energy within the mixing region at a given 
axial station, E m , and the right-hand side is the energy which has entered 
the mixing region up to that station, E z . The two were found to agree 
within 3 percent or better at any point along the axis. This percentage 
was accepted as a sufficiently accurate balance of the energy equation by 
using the assumed velocity profile and temperature-velocity relation. 

4.7 Concentration of Exhaust Oases 

In the development of the governing equations it has been assumed 
that differences in the molecular weights and specific heats of the undi- 
luted exhaust gas and air are negligible. Such differences as may exist 
affect the turbulent mixing process primarily because of their effect on 



1 1 
b * / f V d, i + br i / ? Vi - 

o o 


42 


the density distribution in the constant pressure mixing region. We can 
show that 


P 



(4-66) 


where 

w = W*/Wj 

W* = Molecular weight of mixture at any point 

W* = Molecular weight of undiluted exhaust gas 
3 (at jet exit) 

In general, W will vary through the mixing region if the molecular 
weight of the undiluted exhaust gas and external stream differ. The local 
molecular weight may be expressed in terms of that of the jet, W. and the 
free stream, W TO and of the concentrations of the gases present. If we 
assume that there are effectively only two gases, the jet and the free- 
stream gases, we may represent the mass fraction of the jet gas present at 
any point by c. The quantity W at that point is then 

W - c + (1 - cjW^ (4-67) 


where 

Woo - W>j 

W* « Molecular weight of the free-stream gas 

Specie differences also affect the distribution of energy in the 
mixing region due to associated differences in the specific heats of the 
species. Thus the enthalpy ratios in equation (4-18) cannot be simply 
replaced by corresponding temperature ratios, but must incorporate the 
associated specific heat ratios also. Thus 


where 


P 



p j 
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In exactly the same manner as for the molecular weight, the specific heat 
may be expressed in terms of the local concentration and the specific 
heats of the jet and free stream. 

c « c + (1 - c)c (4-68) 

ir Poo 


where 


c = c * /c* = ratio of c*'s of the free stream 

P oo Poo P n P 

J and the undiluted exhaust gas. 

Consider now the differential equation for the concentration distri- 
bution. In a manner similar to that used for the temperature, we may 
associate the concentration and velocity distributions by the similarity 
of their governing equations. The differential equation governing momentum 
transport in a constant pressure, turbulent axisymmetric flow is: 

Momentum 

p° IS * P v If ■ F & (V u If ) <«-«»> 

The conservation of matter leads to the following equation for the 
concentration 

Species 

P U t! + P V If “ r ch?( r P e m l?) + w (4 " 70) 

where 

w ■ the generation rate of the species associated with c by 
chemical means 

€ m ■ the coefficient of turbulent diffusion for the species 

m involved - eVr U* 

nr n j 

€ ■ the coefficient of turbulent momentum diffusion (eddy 

viscosity coefficient) - e u^ r n U j 

We observe that these two equations for U and c are of identical 
form if w - 0, Hence, we may assume that- 

c - c(U) 
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and obtain from the above equations with w « 0 


(> $ * pv I?)* - 1 [&("*. m ♦ '*•.(£)■ %] 


(4-71) 


Defining the Schmidt number as 


c € 


jm 

u 


and multiplying the momentum equation by dc/dU we may obtain 


A 

dr 


[ r P £ u lS c * 11 Sf] dS + r P £ u S c(l ?) 2 ^ ‘ ° 


(4-72) 


If the Schmidt number is unity, i.e., = e^, we find simply that 


c = A + BU 
c c 


(4-73) 


Hinze and van der Hegge Zijnen (ref . 4-20) investigated experimentally 
the spread of tracer gas introduced into the nozzle fluid of subsonic 
turbulent jets. The turbulent Schmidt number was found to vary slightly 
across these jets and to have an average value of about 1.36 across the 
jet. Such differences compared to the assumption of S c - const. * 1 
should be of little consequence in the present analysis due to its integral 
nature and due to the small difference between the actual Schmidt number 
and unity. 

Using the boundary conditions 

c » 1 
c ■ 0 

the coefficients A c and B c 


when U 


when U ■ U_ 


may be evaluated to yield 


u - u p 
1 - u 


(4-74) 


45 


Introducing the similarity velocity profile 


- * - 
U m 

u - u 

C oo 


where 

U c ■ U c (z) the jet centerline velocity 
n ■ (r - r i )/b 

we obtain 



( 4 - 75 ) 


( 4 - 76 ) 


or 



rj 3 / 2 ) 2 


( 4 - 77 ) 


The mixing calculations provide U c (z), r^z) and b(z) such that 
the entire concentration distribution, c(r,z), maybe determined. 
Calculations for the present situation have been carried out under the 
assumption that Cp^ ■ - 1, that is, neglecting differences between 

specific heat at constant pressure and molecular weight between the free 
stream and the jet. 
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Investigator ( s) 


o 
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o 

s 

a 

o 
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Liepmann and Laufer 
Albertson, et al. 
Corrsin and Uberoi 
pabst 


Warren 

Maydew and Reed 
Baker and Weinstein 
Eggers and Torrence 
Wilder and Hindersinn 


Mach Number Ranges 

= 0 
* 0 

w 0.08 
0.85 

0.69 - 2.6 
0.7 - 1.96 
a 0 

0.47 - 0.94 0.58 

0.94 - 2.8 1.25 
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APPENDIX 4A 


EXTRACTION OF K FROM VARIOUS DATA SOURCES 


In the following the methods are shown by which the eddy viscosity 
coefficient, K, defined by equation (4-48) was obtained from the published 
data of various investigators: Warren, (ref. 4-11) , Corrsin and Uberoi, 

(ref. 4-10) , and Pabst (4-16) : 

Warren investigated the flow field of compressible axisymmetric jets 
exhausting into still air. In his investigation Warren defines the eddy 
viscosity (in present notation) as 


- r *> 


H 

2 


(4A-1) 


where 

r* , * the radius at which the velocity is one-half the 

' centerline (jet) velocity 

We may equate equation (4A-1) with the present definition, equa- 
tion (4-43) since e* must be the same for both. Since U* ■ 0 in 
Warren's study, we see that 


K - l K w 




- r. 


b* 


2 K w\/s 


(4A-2) 


The ratio (r*, g - r£)/b* - r\ x ^ z 
similarity velocity profile, equation (4-25) , to be 


and is found from the assumed 


V* " 0,44 


(4A-3) 


Warren has determined values of based on the definition given 

in equation (4A-1) for his measurements as well as those of Corrsin and 
Uberoi and of Pabst, Equations (4A-2) and (4A-3) were used to calculate 
the K's presented here from those appearing in Warren's report. 

Warren determined his value of by obtaining the best match of 

his analytically selected velocity decay rates to those measured by him. 
In order to correct for differences in Warren's analysis and the present 
one, we also calculated his measured potential core length for his 
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experimental conditions using the present analysis. Another value of K 
was obtained in this manner. It was found that the resultant K's were 
only slightly different from those given by Warren. The calculations are 
tabulated in Table 4A-I. Since these values were more directly comparable 
to the other data obtained, they were used in the correlation. The K's 
computed in this manner for the Corrsin and Uberoi data and for the Pabst 
data were identical to those computed using equation (4A-2) . 

Baker and Weinstein (ref. 4-17) : 

The study of low speed turbulent mixing of co-flowing dissimilar gases 
by Baker and Weinstein also included some measurements for air-air mixing. 
The definition of eddy viscosity U 3 ed in this study was similar to that in 
the present analysis except the mixing zone width, b£ >0s , was defined as 
the distance between points where U ■ 0.05 and U = 0.95. Thus, 


K - ie„b* /b* 

oVV o • 0 5 

" K BW At > 0 .os 

From 

T| - (1 - U / 2 ) 2 /® 

A^q .05 - [ X - (0.05) 1 / 2 ] 2 / 3 - [l - (0.95 ) l/2 ] 2/3 
• 0.844 - 0.086 
- 0.758 

Thus 

K - 0.758 

Mavdew and Reed (ref. 4-7) : 

Maydew and Reed measured velocity profiles and mixing region growth 
in the initial (potential core) region of sub and supersonic jets issuing 
into still air. They compared analytical solutions for the jet velocity 
profiles given by Tollmien (ref. 4-18) , Gortler (ref. 4-19) , and Crane 
(ref, 4-6). These solutions are based on a similarity variable 


£ " oy/Z 
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where 


y is the distance from the point where U = 1/2 

a is the similarity parameter for the velocity profile 
chosen to indicate the spreading rate of the mixing 
region. 


The spreading rate measured by Maydew and Reed is the spread of the 

two points in the velocity profile, where U 2 “ 0.9 and U 2 = 0.1, and is 

denoted as b„ , . 

o . 1 


db 

dz 


A£ 


o» 1 


1 

a 


- At| 


0.1 


db 

dz 


where 


i - S 2 

u ^® 0 . 9 


_ £ 


U 2 =o. l 


- rj 2 
U 2 =0.9 U “0.1 


An c i ■ 

The coefficient K is related to cr and b„ . in the present case by 

O « 1 


K 


(db o /az)‘ 


2o 2 (db/dz) 2(A5„ ,) (db/dz) 

0.1 


K _ AT) o.i db 
2(A^ 0>1 ) 2 dz 

For the Crane velocity profile used by Maydew and Reed 


(4A-4) 


A| 0 


. i 


1.3 


For the present velocity profile 

At ) 0<1 " 0,333 


Thus, equation (4A-4) becomes 
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(4A-5) 


K - 0.0905 (db /dz) 

0*1 

which was used to obtain K from these data-* 

Eggers and Torrence (ref. 4-13) : 

This study is concerned with the mixing of an axisymmetric jet with 
a co-flowing stream. The jet Mach- number ranges from 0.47 to 0.94 and the 
f ree .-stream Mach number ranges from 0.58 to 1.3. Both streams are air and 
have the same stagnation temperature. 

The eddy viscosity definition used in this investigation was not based 
on the velocity differenceJbetween_the_sbrfiams_.but only on the jet velocity. 
Thus, 

e* - K_b* U* (4A-6) 

TS 0.S5 c 


where 

b - r* (U - 0.5) - r*(U - 0.95) (4A-7) 

0.95 

Comparing equations (4A-6) and (4A-7) with equations (4-49) and (4-46) 
respectively, it is found that 

K ■ - "-’I 

and 

- °- j55 


Thus, 

K - 0.355 - Uj (4A-8) 

Equation (4A-8) was used to calculate K from the K £ of Eggers and 
Torrence. 

Wilder and Hinderslnn (ref. 4-1) : 

In this investigation the turbulent mixing of two-dimensional com- 
pressible air streams was comprehensively investigated. Stagnation 
temperature ratios, To w , of 0.5, 0.616 and 0.803 were studied within the 
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Mach number ranges 0,94 <_ Mj 2.8 and 1.25 <, 1.6. Stagnation 

temperature profiles were obtained for a total of 19 conditions within 
these ranges of Mach number and stagnation temperature ratio. The eddy 
viscosity was evaluated for each of these conditions based on an error 
function ^velocity profile. The Crocco relation, equation ( 4 — 20 ) , was 
assumed to hold between the stagnation temperature and the velocity fields. 
Therefore, the error function profile was also assumed to represent the 
stagnation temperature profile, Thus, it was assumed that 

B = ~ (1 + 0(C)) (4A-9) 


where 


0(0 



(4A-10) 


C - (rj. 5 - r*)/y(2€ o /U*L) z* (4A-11) 

B “ (T* - T* ) / (T* - T* ) - (T - T )/(l “ T ) (4A-12) 

° °eo °j °oo ° °eo °oo 

r* ■= the radius at which B = 0.5 (4A-13) 

0.5 

z* « z + a (4A-14) 


a = a displacement of the origin to account for heat transfer 
between the streams prior to the jet nozzle exit. 

The eddy viscosity was formulated simply as 

e* - f (4A-15) 

After "a" was selected to provide axial similarity for the profiles, 
values of (6 Q /L) were determined for 19 conditions by seeking the best 
fit of the data to equation (4A-9) . The selection of "a" appeared to be 
somewhat arbitrary so it was decided to eliminate z in computing K 
from (Cq/L). The following is the calculation procedure which appeared 
most satisfactory although three others were used and generally yielded 
results within 10 to 15 percent of these. 
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By comparing equations (4A-15) and (4-48) , it is seen that 


e o z U q° 

U*L b ^ “ U J 

00 


(4A-16) 


We also note that matching the slope of 0 at its maximum (i.e. , at 
0-0.5, £ - 0) we obtain a good match of the 0 profile with the theory. 
It may be shown from equations (4A-9) and (4A-10) that 




Y IT 


(4A-17) 


From equations (4A-11) and (4-27) it is seen that 


ir o z 


(4A-18) 


where ti - t) at 0-0.5, from which we obtain 


& _ _ 


!o m 


U*L b dT l 
00 


(4A-19) 


Using (4A-17) and (4A-19) it may be shown that 



. <&! 


(4A-20) 


Equations (4A-16) and (4A-20) are now combined to yield 


2ttU*L 

00 


_ 

" dtj 


(4A-21) 


After evaluating (d0/dt] ) we may use the data from Wilder and 

T)o • 5 

Hindersinn directly in equation (4A-21) to determine K. We now proceed 
to evaluate (d0/d^) using the relation 
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dg. _ d£ dU 
dT l dU dT > 

j?rom equation (4-23) for U c = 1 (all data are for near field conditions) . 

r (1 " U J 2 I- (1 - U J _ 

0 - 1 - D C2 -y - _ - T I U + DC® y _' T if® (4A-22) 

^ °oo °oc 


where 


c j " < u j S / 2 >/ H j “ M j /C 1 + x t Jl M j ) 


-c 


0 for Crocco relation 

for modified Crocco relation 


Prom 


dg. 

dU 


(1 - UJ : 


* 1 + D C® i _ T (2U - 1) 


(4A-23) 


From equation (4-25) 


we observe that 


dU 

d^ 


U - (1 - n 3/2 ) 2 


- 3 (1 - t, 3/2 ) 1/2 


(4A-25) 


or, using equation (4A-25) 


- Ji'/’d - 


(4A-2S) 


Thus, combining equations (4A-24) and (4A-25) 




^ - - 3 u 1/a (l - u 1/2 ) l/3 |i + d ^ (2U - 1)| 


(4A-26) 


We see that for the Crocco relation 0 - U and (d0/drj) - (du/drj) 


dg. 

d n 


dg. 


du 

dr\ 


'o.s 


%. 5 


'0.5 


- 1.41 


- 3(0.5) l/a (l - 0.5 l/2 ) l/3 


for Crocco relation 


(4A-27) 


For the modified Crocco relation, equation (4A-22) may be solved for 

1. 


U with D 

o 5 


U 


0.5 


1 - 


/ / 1 - T °, X 
- «»>* + v 1+ \ o j (i - u jv 


(4A-2B) 


Solution of equation (4A-28) for the particular values of Mj, T Og0 » 
and U, and substitution into equation (4A-26) yield (de/dt^ 5 ) to be 
used in computing K from equation (4A-21) . The basic data and tabulated 
computations are shown in Table 4A-II. 

A tabulation of all final data appearing in figure 2 may be found in 
Table 4A-I and Table 4A-II. 
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Experiment 


0.69 1.0 10.2 
0.69 0.7 9.4 
0.97 1.0 12.0 
0.97 0.7 10.7 
1.51 1.0 14.1 
2.60 1.0 25.0 
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5.0 INITIAL INTERACTION BETWEEN TRAILING 
VORTICES AND JETS 


5.1 Description of Model 

During the initial stage of the jet exhaust-trailing vortex 
interaction process, the hot gas jets retain their integrity and 
axisymmetric shape while the trailing vortex sheet wraps up into a pair 
of concentrated vortices. Because of aircraft angle of attack, the 
jets originally cant downward and tend to move away from the vortices 
in the crossflow plane. The component of free-stream velocity in the 
crossflow plane causes an upwash about the jets which, due to crossflow 
drag, turns the jets into the free-stream direction. The downwash on 
the jets caused by the trailing vortex sheet diminishes this effect 
somewhat. 

In determining this initial interaction between the trailing 
vortices and the jet, it is assumed that the jet acts as a solid boundary 
expanding outward at a rate given by the characteristics of the turbulent 
jet alone and that the crossflow drag coefficient is uniform and given 
in accordance with the work of Wooler (ref. 5-1) . It is also assumed 
that the jet remains axisymmetric about its curving centerline, that its 
other properties are the same as for no curvature, and that the jet boundary 
is the proper boundary to image the external vortex sheet. Both sides 
of the wing have been accounted for in the rolling-up calculation. 

In order to illustrate the method, consider the two exhaust jets 
together with one pair of external vortices, as shown below. 



Vertical plane of symmetry 




The image vortices inside each jet are included. We note that the 
sum of these internal image vortices is zero for each jet, so that the 
influence of the left jet on the right vortex and right jet can be 
neglected. However, the effect of the left vortex and its image in the 
right cylinder on the motion of the right vortex is of higher order and 
must be included. A total of 10 equal strength vortices on each wing 
panel, located to correspond to an elliptical loading at the trailing 
edge, have been taken to model the vortex sheet, which is constrained to 
have mirror symmetry. In addition to the pressure forces on the jet due 
to the vortices, and the crossflow drag force due to the initial down- 
ward cant of the jet, there will be a net buoyant force on the jet 
due to the temperature difference between the hot exhaust gases and the 
ambient air. The initial interaction of the jet and trailing vortices 
and subsequent motion is the result of all of these combined effects, 
and the governing differential equations for the trajectories are derived 
in the following section. 


5.2 Differential Equations for Jet-Trailing 
Vortex Trajectories 

The differential equations for the mutually interactive motion of 
the jet and trailing vortex sheet are derived by considering the momentum 
equations for a thin cross-sectional slice of the jet exhaust under the 
combined action of the pressure forces due to the vortex sheet, crossflow 
drag, and buoyancy. Consider the side and top views of the right-wing 
vortex centroid and jet as sketched below 



z 
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The jet is initially canted downward as viewed in the yz plane at the 
thrust axis inclination. As the jet proceeds downstream, it ingests 
free-stream air at a rate m, given by the turbulent jet characteristics. 

a 

This tends to align the jet with the free stream, as does the action of 
the crossflow drag force per unit length, F^. The action of the net 
buoyant force Fg per unit length tends to cause the jet to rise verti- 
cally, while the influence of the trailing vortex sheet pressure forces 
(p^ J p ) is primarily such as to cause the jet to move in and around 
the center of vorticity. 

Thus, the nondimens ional equations of motion of the thin slice of 
the jet are given by 



(5-1) 

(5-2) 


5.2.1 Vortex forces .- The vortex forces per unit length (F^.F^) 
can be obtained by integrating the pressure in the crossflow plane, as 
given by the unsteady Bernoulli equation, around the boundary of the 
moving jet. 

This unsteady surface pressure can be calculated as follows. Consider 
the crossflow plane view of the moving cylinder of exhaust gas and a 
typical pair of vortices centered with respect to the plane of symmetry 
of the aircraft. 



The complex velocity potential for this external vortex pair plus 
images with respect to the moving cylinder at ( x c »y c ) cylinder 
coordinates Z 1 « + iy* is given by 


W(Z X ) 
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(5-3) 


where V r is the radial velocity at the boundary of the jet as given 
by turbulent jet theory. 

This is easily generalized to N vortex pairs by 


W(Z X ) 
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(x c + iy^) 


r e v r m Z X 


(5-4) 


The pressure is determined from the unsteady Bernoulli equation, which 
is given in nondimensional form in the accelerating (x^y^ coordinate 
system by Lamb (ref # 5-2, p. 233) as 


Poo “ Pc 


\ (v 2 
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(5-5) 
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where time derivatives have been replaced by their equivalent longitudinal 
derivatives, i.e. <>/dt - U w d/dz and 0, v, and w are determined 
from W(Z X ) in the usual way, i.e. 

0 - R.P. W(Z X ) (5-6) 


and 


dw 

dZ, 


v - iw 


(5-7) 


Finally, the forces (F^F ) ou the cylinder due to the vortical 
flow field described by W(Z X ) are determined by 

2tt 


vx 


+ iP vy" - / r eP <r e )e10 


de 


(5-8) 


where p(r ) is the nondimens ion al pressure at the jet boundary 
e 

r - r g as given by equation (5-5) . 

While the expression for the complex potential (5-4) is compact, 
the operations indicated by equations (5-5) through (5-8) are rather 
lengthy. The resulting expressions for ^Fy^Fyy) are given by 



(continued on next page) 
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5,2.2 Crossflow drag force .- The nondimensional crossflow drag 
force per unit length of the jet is given by the relation 


F_ % 
Dx x 


F % 

Dy y 


c s i p- 1^1 


(5-13) 


where C D is the crossflow drag coefficient per unit length and is 
given from the work of Wooler (ref. 5-1) by ■ 1.8, and xfj, is the 
total nondimensional crossflow velocity (free stream plus vortex induced) 
normal to the jet centerline. For smal-1 jet inclinations to the free- 
st re am direction, is given by 


1 *( XI v vx - u « ijtK + ( L 
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all external 
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vy 




(5-14) 

where the summations Taj and Hv indicate the total external 

vx vy 

vortex-induced velocities at the jet centerline in the x and y 
directions, respectively. For example, for a single vortex pair 
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(5-15) 


(5-16) 


5.2.3 Buoyant force .- The buoyant force per unit length Fg 
represents the net vertical force on the jet due to the temperature 
differential between- the jet and the surrounding ambient air, and is 
equal to the difference in the upward force due to the atmospheric 
pressure gradient acting on the cylinder and the downward force due to 
gravitational attraction. In nondimensional form, Fg is given by 


f B " vr e { P* 




★ 2 
j 


(5-17) 
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5.2.4 Equations for iet and concentrated vortex pair .- Although 
equations (5-1) through (5-17) are somewhat lengthy for the general case 
of N external vortex pairs, they simplify into more reasonable forms 
for the case where the vortex sheet can be adequately represented by a 
pair of concentrated vortices. The equations of motion for the jet 
centerline and the vortex centroid trajectories are 
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and where | ^ | is found from equations (5-14) to (5-16). 
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5.3 Initial Interaction Calculations 

The equations of motion given by equations (5-1) through (5-17) 
have been integrated with 10 equal strength vortices distributed on 
the trailing edge of each wing panel according to an elliptical loading, 
and for the initial conditions and physical parameters corresponding to 
the calculative example discussed in- Section 9. Comparisons were made 
in the motions of the cylinder and vortex centroids for this N ■ 10 
case with the analogous results- obtained by concentrating the vorticity 
into a pair of vortices located at- the center-of-gravity of vorticity 
(i.e. integrating eqs. (5-18) through (5-21)), and~-the results were 
virtually identical. Consequently, all the calculations of the initial 
jet-vortex interaction stage and also the ingestion stage were performed 
with the simpler concentrated vortex-pair model. 

The interaction calculations were carried out according to the 
following criteria. The equations of motion were integrated from the 
aircraft trailing edge downstream to a point where the outer jet boundary 
(which is considered as solid during this phase) has expanded to a size 
sufficiently large that the image vortices inside the jet begin to 
interact with the near external vortex centroid in a strong and spurious 
manner. At that location, the initial interaction calculation is termi- 
nated since- it can now be safely assumed that the external vortex has 
begun to tear apart and ingest at least the closest portion of the hot 
jet. The jet disintegration and ingestion process, described in the 
following section, is then initiated, using as starting conditions, the 
locations and crossflow velocities predicted by the initial interaction 
model. 


6.0 INGESTION OP JETS BY TRAILING VORTICES 


6.1 Description of Model 

When the jet and trailing vortex begin to interact strongly, the 
model constructed for the initial phase of the interaction process is 
no longer valid and must be modified. In particular, the assumption 
of the outer jet boundary as being solid and impervious to the crossflow 
flow field must be relaxed in order to track the various jet segments as 
the jet disintegrates and is ingested into the vortex recirculation 

region. 

To model the break-up and ingestion process, the technique used is 
to track the entire jet using the initial interaction equations of 
Section 5 to the point of strong jet-vortex interaction, then subdivide 
the jet into a number of equal area segments (5 to 10) , initialize the 
various segments with crossflow velocities predicted by the initial jet- 
vortex interaction model, and subsequently track the various jet segments 
under the combined action of the trail ing-vortex pair pressure forces 
plus buoyancy effects. During this process, it is assumed that the jet 
segments are sufficiently small that the mutual interference between 
segments is negligible, and also that the segments can be considered 
to have no internal structure (i.e., do not satisfy a no-flow boundary 
condition at their surfaces) so that there are no image vortices nor 
dipoles located inside them. In addition, since the jet is now 
essentially aligned longitudinally with the trailing-vortex pair, the 
crossflow drag can be neglected. 

6.2 Differential Equations for Jet Disintegration 

The differential equations governing the various segments of the 
jet during the break-up process can be obtained either from first 
principles, or as a limiting case of the initial interaction equations 
as both the radius and rate of change of the radius of the segments goes 
to zero. Also, because the mutaal interaction between the segments is 
assumed to be small, the motion of each individual segment can be calcu- 
lated independently of the others. Consequently, the equations governing 
the location of an arbitrary segment of the disintegrating jet are given by 
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U «A dz + dz J 


(6-4) 


Equations (6-3) and (6-4) indicate that the motion of the concen- 
trated vortex pair is strictly vertical at a downward rate equal to 
their mutually-induced velocity. This, of course, is expected since 
the assumption was made that -the individual jet segments are sufficiently 
small so that their influence on the motion of the trailing-vortex 
pair is negligible. 


6.3 Disintegration Calculations 

The jet disintegration trajectory equations (6-1) through (6-4) 
were integrated according to the following criteria. At the end of the 
initial interaction phase, the jet is divided into a number of equal 
area segments (five were found to be sufficient) . The locations and 
crossflow velocities of the jet segments predicted by the initial inter- 
action calculated are used as starting conditions for equations (6-1) 
through (6-4) . These equations are then integrated downstream for each 
segment of the jet to a point sufficiently far that it is clear that the 
particular segment is either captured or escapes. The criteria used 
for capture is that the trajectory of the segment becomes either a 
closed curve about the vortex or a spiral into the core. 

To investigate the effect of various initial velocities of the jet 
segments, a series of disintegration calculations were carried out 
beginning at various longitudinal positions upstream of the-point where 
the initial interaction calculations are terminated. 

It was found that for the typical parameters and conditions corres- 
ponding to the calculative example, the vortex strength is sufficiently 
high to capture all of the hot jet exhaust within a relatively short 
distance downstream. This result was relatively insensitive to the 
starting position of the disintegrating phase. These calculations are 
described in detail in Section 9. 
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7.0 EQUILIBRIUM VORTEX STAGE 


7.1 Introduction 

It has been calculated that all the exhaust products from the engines 
are entrained in the recirculating flow. The actual distribution of 
temperature in the recirculating flow is not necessarily an indication of 
the engine exhaust gas concentration as during the jet phase because the 
vortices cause radial pressure gradients which themselves cause temperature 
differences even when no exhaust gas products are present. It appears that 
the turbulent transport processes in the recirculation region will control 
the distribution of temperature within that region if the stage lasts 
sufficiently long. 

Within the recirculation region the field near the vortices will be 
dominated by the local vortex and will take on the character. of an equili- 
brium turbulent vortex. In the outer reaches of the recirculation region 
between the vortices, the effect of both vortices will be present. It is 
possible to obtain a good idea of the characteristics of much of the flow 
by considering it to be simply two isolated equilibrium turbulent vortices. 
This is particularly true of the core region where most of the exhaust gases 
tend to be entrained. Accordingly we will now consider the questions of 
the temperature and concentration distributions for an isolated equilibrium 
vortex. 

7.2 Kinematics of an Equilibrium Turbulent Vortex 

Before determining the temperature distribution of an equilibrium 
turbulent vortex, if is helpful to describe the regions of the vortex and 
its tangential velocity and circulation distributions. A turbulent vortex 
can be considered to fall into three parts as discussed by Nielsen and 
Schwind in refereice 7-1. The center of a turbulent vortex is in solid 
body rotation with the tangential velocity proportional to the radius. 

This region merges into the central part of the vortex wherein the circu- 
lation varies logarithmically with the radius. In this region the tangen- 
tial velocity reaches maximum v x at a radius x x . Beyond the logarithmic 
regi«-> of the vortex is a "wake" region analogous to the weJce region of a 
turbulent boundary layer. In this region the circulation decays exponen- 
tially with radius. 
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With this brief description, we now introduce figure 7-1 taken from 
reference 7-1. In this figure x L is the radius of the "eye" of the 
vortex in solid body rotation, r.. is the radius of the outer edge of the 
logarithmic region, and r Q is the outer edge of the vortex taken where 
the circulation is 99 percent of the total circulation. The following 
empirical results taken from reference 7-1 are based on the data of 
Hoffman and Joubert (ref. 7-=2) and give the circulation distributions in 

the three regions. 


Eye: 

H'i'i 

II 

*1* 


(7-1) 

Logarithmic, 

Region 

Ci < r < 

K 

Ki 

0.938 In (j^) * 1 

(7-2) 

Wake: 

i r S- r o 
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(7-4) 
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(7-5) 


r-L = core radius based on experiment 


v x * tangential velocity at core radius 


* v i r i 

The tangential velocity distribution is simply given by v - K/r. 

7.3 Vortex Temperature Distribution 
We now determine the temperature distribution through an equilibrium 
turbulent vortex. 

7.3.1 Vortex "eve".- Consider first the "eye" of the vortex in solid 
body rotation. For this region it is known that rotation has a strong 
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influence in damping the turbulent fluctuations so that to the first 
approximation we will treat this region as laminar. 

The governing equations as taken from Hall (ref. 7-3) are 
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where the total enthalpy is 
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I ■ c T + i (u 2 + v 2 + w 2 ) 
o p ^ 


a ■ Prandtl number 


Assumptions in deriving these equations include 

(a) Viscosity and thermal conductivity are constant 

(b) Laminar flow 

(c) Axial symmetry 

(d) Axial time derivatives are small compared to radial time 

derivatives 
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We now add the following assumptions 


(e) o-l 

(f) Steady flow 

so that .equation (7-10) becomes 
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(7-11) 


Trailing vortices in the equilibrium vortex stage are many miles long 
and their core diameter changes very slowly. Thus u and Sl Q /Sz are 
very small. We will take them to be zero so that 

“ M-v 2 j = 0 (7-12) 

This equation simply states that the work done on the eye of the vortex 
by viscous shear is transferred back out as heat. 

We now have 
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Since 

v - 0 




C x - 0 


and 




ft 2 r 2 

2 


+ C 2 


Since u « v, w « v 


(7-14) 


c p T + 


v£ 

2 


n g r 2 

2 


+ Cg 


79 


r 


8 



Let be the temperature at the outer edge of the edge of the eye of 

the vortex at radius r^. 


r ■ r. T ■ T. v ■ v. • fir. 

i i ii 

(7-15) 

f! ' 

ii 

The resulting equation is 
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T ■ T i r < r i 

(7-16) 

h 

For a Prandtl number different from unity, equation (7—12) yields 
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(7-17) 


We will use the simple a ■ 1 solution with the eye of the vortex 
solid body rotation at uniform static temperature. 
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(a) Steady flow 
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(c) Density fluctuations are zero 

(d) Axial symmetry 
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The resulting equations deduced fror,t reference 7-3 are 
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At this point we introduce the Reynolds stress due to turbulence 

and the eor r esponding heat transfer 
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In rotary flow the shearing stress is taken proportional to the eddy 
viscosity as follows 
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with the kinematic term zero for solid body rotation. With respect to 
the eddy thermal conductivity, it has been suggested (ref. 7-4) that 
adiabatic temperature changes associated with radial pressure gradients 
do not act as a driving_force for heal, transfer 30 that 



Now for adiabatic changes 
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For small fluctuations we hc.ve 


I' ■ c T' + vv' 
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(7-24) 
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(7-27) 


u'l' ■ c u*T' + v u'v' 
O p 


incorporating this result into equation (7-21) and letting u - 0 yield 
£ <rp - £ [rp(c p *.v ^>] - ^ I? 

- je + v v ^ ■ ° 

2 S? V' pr- t r J J 

tr 

or assuming the turbulent Prandtl number is unity 

■h ( rp v t [ c p & 


(7-28).- 
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or 


& (V + T-) - 
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(7-30) 


The right-hand side of the equation can be treated by„_use of 
equations (7-20) and (7-23) 

& <**p - - & [^ pv t(lr - Oil 0 


(7-31) 


Thus 


Tpvp - °* r (t^ ■ 0 

substituting into equation (7*"30) , we find 

& (v + f) * ^ C = C * r (£ - r) 


(7-32) 


(7-33) 
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In seeking solutions to equation (7-33) we need to specify the radial 
distribution of the tangential velocity, v(r). Several possibilities 
occur. We can specify the radial distribution for an equilibrium vortex 
as given by equations (7-1) to (7-3) . For a simple solution we might 
specify a variation corresponding to a potential vortex in the field 
external to the eye of the vortex. It is of interest to compare equa- 
tion (7-33) with equation (7-14) before proceeding with any solutions. 

The first thing to- note is that (dv/dr - v/r) is zero for solid body 
rotation so that this term in equation (7-33) does not contribute to 
dT/dr at r - r jL . It is also noted that a factor of two appears before 
the v 2 /r term in equation (7-33) while a factor of unity appears in 
equation (7-14) . There is thus a discontinuity in the temperature gradient 
at r ■ r^ due to the different heat transfer laws for laminar and turbu- 
lent flow. To alleviate this difficulty requires a buffer region between 
the laminar and turbulent region which has not been introduced into the 
theory. However, this simplification of the theory should not interfere 
with obtaining some interesting approximate solutions. 

For the assumption of a potential vortex external to the eye of the 
vortex, a particularly simple solution can be obtained. In this case the 
tangential velocity is given by 


v 



(7-34) 


The term proportional to C 3 C 4 in equation (7-J3) must be zero since it 
will produce a log r term which is not regular a3 r — Thus we obtain 
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If we take as a boundary condition 

T - T* 
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then 
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(7-37) 


The temperature of the eye is then 


The eye will be colder than the free stream. This result should give an 
approximate idea of how cold— the eye will be for the case of-an equilibrium 
turbulent vortex. 

7.3.3 Improved solution based on logarit hmic circulation distri- 
bution.- For arv equilibrium vortex the circulation distribution given by 
the logarithmic law yields a good approximation even in the wake region 
even though it does not have the desired properties as r — ». If we 
limit the solution to the region r ^ r^, we can use the logarithmic 
circulation distribution to obtain an improved solution. We can retain 
both terms in the right-hand side of equation (7-33) since the range of 
r is limited. At r - r Q we can apply both temperature and temperature 
gradient boundary conditions assuming that the stagnation enthalpy is 
constant for r > r Q 


&(v + *f) m 0: 

Subject to this outer condition, it is readily found from equation (7-33) 
that 

c T + ^ - c T* - 2 {(0.928 K 1 )2[F 1 (r <J ) -F,(r)J 

r 4 r 

+ 0.928 - 0.928 In ri)[F 2 (r 0 ) - F z (r)] 

+ Ki(l - 0.928 In r x ) 13 |f 3 (r Q ) - F 3 (r)]} 

+ C s [o.928 K x (In r - in r Q ) - 0.928 K x (ln*r - ln a r Q ^ 

- 2K X (1 - 0.928 In ri) (In r - In r D )] (7-40) 


r > r Q (7-38) 

r > r G (7-39) 


given by 



84 




WJffl 


Tiir»r 




8.0 ASYMPTOTIC BUOYANCY PHASE 


8.1 introduction 

It is of interest to know the ultimate disposition of the exhaust 
gas products in the stratosphere. Several studies of this question are 
available, that of Overcamp and Pay (ref. 8-1) . In reference 8-1 smoke- 
stack plume theory is applied to the hot wake of an aircraft. In ref- 
erence 8-2, several dispersive mechanisms are discussed; namely, wake 
collapse, wind shear, buoyancy rise, and turbulent- diffusion. Without 
discussing these results in detail, we note that they are not mutually 
exclusive, not linearly superposable, and certain ones are dependent on 
atmospheric turbulence and winds. There is noticeable difference between 
the buoyancy results of the two groups of authors. 

In this section are set the tasks of finding the asymptotic solution 
for a quiescient atmosphere. We will address the problem of a hot region 
of recirculating flow associated with the pair of turbulent vortices. We 
will assume that at some distance behind the aircraft the axial flow in 
the recirculating region had died out, that the vortices are close to 
horizontal, and that the wake moves along with any horizontal wind. 

The paper of Fay and Overcamp concerning the effects of- buoyancy ©n 
the width of the recirculation region includes a correlation of the width 
of such regions obtained- from contrail data. The paper also includes an 
interesting comparison between the data correlation and an approximate 
solution obtained from the stack plume theory of Hewett, Fay, and Hoult 
(ref. 8-4) . In order to assess the theoretical basis of this comparison 
we now examine the theory of that reference. 

Consider a hot plume emerging from a stack in a uniform cross wind 
so that the plume is stationary in stack coordinates. Let v be the 
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horizontal wind velocity and u the uniform axial velocity of the gas in 
the plume at angle 0 to the horizontal. The rate of entrainment € is 
taken proportional to the relative wind velocities parallel and perpen- 
dicular to the plume axis as follows 

e - 2irR[a|u - v cos d\ + 0 |v sin $|] (8-1) 

In this equation a and 8 are constants of the order of unity. Subject 
to the condition that R ■ 0 at x • 0 (at the stack) the following 
approximate solution is obtained from reference 8—4 for the plume radius 
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O) - 9 - Brunt-Vaisalla frequency (8-3) 
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Expansion of the above expression for small values of x/i^S yields 
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Consider a distance d^ 
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such that the cumulative buoyancy force 
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and introduce a new reference length l such that jj 

l - UbdJlV* (8-5) 

Then we obtain 

Of) ■ ifYXf^T (a - 6 > II 
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This power law is compared with the contrail data of reference 8—1 in fig- 
ure. 8-1 assuming 8-1. After an initial period, the contrail wake grows 
in size at a rate equal to or greater than that given by the above power 
law. 

Examination of the limitations on x in equation (8-2) shows that 
the power law is valid only for very short times. We can write 



cat < ir 


For T = 392°R we find that 


cu - 0.021 rad/sec ; 

and the maximum time corresponding to tt radius is about 150 seconds. 

Since the power law is valid only for much lesser times than the maximum 
time, it should be valid only for periods less than a minute. What is of 
interest is that the power law is in fair agreement with the contrail data 
for much longer periods of time. However, we cannot conclude that buoyancy 
is the dominant mechanism since the buoyancy theory has been extended far 
beyond its range of validity. This observation is not necessarily true 
for a supersonic airplane. 

In this section we carry out a detailed derivation using the above 
entrainment model for a hot wake which is initially horizontal and which 
possesses no axial nor vertical velocity. Kinetic energy effects are 
initially included, and the case of stratification associated with a 
uniform atmospheric temperature is treated. An analytic solution in 
closed form is found. The complete mathematical-solution is not carried 
out, but numerical results obtained on a computer are presented for an 
initial value of R not equal to zero. These are compared with the data 
correlations and theory of Fay and Overcamp. 

8.2 Horizontal Wake Analysis 

It is to be expected that the axial flow in the recirculation region 
will damp out some distance behind the aircraft, but that the vorticity 
and heat cohtent of the recirculating flow will still be substantially 
unchanged. We assume that tne wake is nearly horizontal with small vertical 
velocity at some point in its equilibrium vortex phase. Consider a crons 
section of this wake at time t - t^ 
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Let the average density in the wake be p. Let h be distance measured 
vertically from the t - t L altitude. The atmospheric density is p*, 
and at t - t £ is designated p^. Let the vertical rate of rise of the 
wake be w. We will assume the wake is circular of radius R for simpli 
city. Let m i be the mass of air in a section of wake of thickness Ax 

at t = t^ 

m^ - ttR| p^ Ax (8-7) 

Let m g be the time rate of mass entrainment into the wake section due 
to its rate of rise w. Thus we have 

m g ■ 2vR Ax alwjp^ (8-8) 

where a is a dimensionless entrainment coefficient. The total mass in 
the wake section at time t is then 


m 


m^ + 2 tt Ax a 


/ 

t i 


Rw Poo dt 


(8-9) 


Since we will be dealing only with positive values of w,__W£_ remove the 
magnitude sign. 

The vertical momentum equation can be obtained by equating the 
buoyant force equal to the rate of change of vertical momentum. The 
buoyant force P fi is 

Fg ■ 7 rR s 'Ax(p 06 - p*) g (8-10) 


With 
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we find that 
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dt 
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( 8 - 12 ) 


with the initial condition 


dw 

dt 


3 ( Poo. “ Pi) “ 2 “ W lpco/ R i 

Pi 


(8-13) 


The application of the law of the conservation of energy to the 
problem is of particular interest because it yieldb a closed solution 
for the size of the wake cross section as a function of altitude. The 
change in total enthalpy of the recirculation fluid between time "t" 
and "t^ 11 is equal to the sum of the work done on the wake by the buoyant 
force during this time, the energy added to it by the entrained mass, and 
the work done by the fluid against atmospheric pressure. 
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( c p ? + t) ■ m i (v'i + t) ■ / 1 


F b w dt + m e c p T 00 + 


f i 


dp 


(8-14) 


In this equation it is not necessary to include potential energy terms 
if we use the net buoyant force in calculating the work done on the wake. 
This is true because the increase in potential energy between two points 
a short distance apart, d (mgh) , is equal to work done on the wake (rng)dh 
plus potential energy (gh)dm added to wake. 

Here we assume that the external air temperature is constant with 
altitude at T*. The temperature T is the mass averaged temperature of 
the recirculating flow and is given by 


T - 
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( 8 - 16 ) 


where R is the gas constant for the recirculating flow. Since 
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we have 
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In differential form equation (8-14) can be written 
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Neglecting w 2 compared to c T* anc ^ assuming w is positive yields 


(2a C p T* + - 2CpT*dR 
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The solution of the equation subject to the initial condition R * R^ 
when h ■ 0 is 


a + 


9R 


2y*g< 


G 

> 2aR_T* / 


gRi N bg/2YR g T 


(e 


V 

Fo^; small values of hg/2YR g T* , the approximate solution is 
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R ■ R i + uh 


(8-23) 


We have assumed in the derivation that the average pressure of the recir- 
culating flow is equal to the atmospheric pressure at the height of its 


center. 


8.3 Solution of Vertical Momentum Equation 


We can write equation (8-12) in the form 
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(8-25) 
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We will determine the explicit functions F(h) and G(h) and thus show 

that equation (8-25) Ls a linear equation for w* . 

Starting with the integral, we have 


t n 

j Rwp,o dt ■ J R P w dz 


(8-28) 


For an isothermal atmosphere of temperature T* > we have 
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Carrying out the integration yields 
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(8-35) 


A formal mathematical solution can be obtained in the form 

. e-/ F(h)dh [/G( f| )(e/ F(h)dh )d,, ♦ c] 

subject to the initial condition 

w ■ 0 at t * ti 

We have chosen to do the solution numerically for our case. 

8.4 Calculative Example for a Supersonic Airplane 
and Comparisons with Data and Theory 

Once the wake is a short distance behind a supersonic airplane, it 
has subsonic velocity with respect to an observer on the ground. It is 
thus expected that a subsonic buoyant wake theory should be applicable to 
the supersonic airplane case. For the purpose of illustration, we now 
present two example calculations for an example supersonic airplane using 
the following parameters : 

Initial altitude: 70,000 ft (21.3 km.) 

Flight speed: 2,670 ft/sec (815 m./sec) 

Wake temperature rise: 9.3°F (5.16°C) 

Free-stream temperature: T* - 392. 3°R (217. 9°K) 

Entrainment coefficients: a ■ 0.3 and a ■ 0.9 

Initial radius: Ri ■ 42.7 ft (13.0 m.) 

Initial density: p W;L - 0.0001389 sluq/ft 3 

Ratio of specific heats: y * 1.4 

The various characteristic lengths used in the Fay-Overcamp theory 
of reference 8-1 are estimated to be 

£ b - 1.90X10 -4 ft 
■» 163,500 ft 
£ - 172.5 ft 

The calculations have been made by determining R as a function of 
h from equation (8-22) . The vertical velocity w was determined as a 
function of h by integrating equation (8-25) numerically. The time 
scale has been established by the equation 
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The radius R has been nondimensionlized in the form 2R /i and the 
distance x in che form x/d b so that the results can be presented 
on figure 8-1 for comparison on the -d at a and the Fay-Overcamp theory. 

First the comparison of the present calculations with the theory of 
Fay. and Overcamp, equation (8-6), shows generally fair agreement within the 
middle~of- the range. The slope of the Fay and -Overcamp theory is slightly 
larger in the mid-range- At small values of x/d b there is a divergence 
between the theories because of the difference in initial conditions. For 
large values of x/d b the present calculations show a peak at approximately 
x/d b --2.5. The Brunt-Vaisalla frequency for the present case of 0.021 
radians/second yields a value of x/c^ - 2.45 for v radius. The theo- 
retical solution thus peaks at just about this value. However} the power 
law, equation (8-6) , does not exhibit a peak value of R. It is noted 
that the final radius of the-wake is about 50 percent greater for a - 0.9 — 
than for a - 0.3. 

With- regard to comparison between the present calculations and the 
correlated data, several points are significant. In the intermediate 
range the slope of the data points is about twice that of the theory. 

One point of interest is that the values of x/d^ corresponding to ir 
radians for the data of figure 8-1 varies between 0.62 to 0.84. Yet no 
maximum value of R occurs experimentally for values of x/d^ much 
greater than these values. It thus appears that there are significant 
physical factors not accounted for in the present buoyancy model. It is 
also of interest that- the- value of x/d^ for ir radius is 3 to 4 times 
larger for the supersonic case considered here than for the subsonic 
cases shown. 

The results for a - 0.3 and a - 0.9 show about a four-fold and a 
six- fold overall increase in wake radius at the peaks, respectively, 
due to wake entrainment. The value of o - 0.3 corresponds to a com- 
pletely quiescient atmosphere. This value was obtained from -the work of 
Tsang on line thermals described in reference 8-5.— By placing a cylinder 
containing heavier fluid horizontally in a tank of lighter fluid and 
suddenly withdrawing the cylinder horizontally he created an initial 
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condition of a cylindrical mass of fluid of higher density in a surrounding 
fluid of lower density without disturbing the fluid. His mean value of a 
obtained from his figure 10 yields « - 0.33. Richards, reference 6-6, 
performed similar tests but did not take sufficient care to minimize 
initial disturbances. His values of <i range from 0.316 to 0.67. It 
is thus apparent that the value of a can be doubled by moderate initial 
disturbances. We have used a * 0.9 as an approximation for the real 
stratosphere . 

There is evidence, as pointed out above, of serious shortcomings in 
present buoyant wake models. One factor which has received scant attention 
is the role of changing vortex strength during the buoyant rise. During 
entrainment circulation is induced in the recirculation region of opposite 
sign to that existing just behind the aircraft. Account should be taken 
of this phenomenon in the fluid mechanical model and possibly in the energy 
equation. 

Overcamp and Fay (ref. 8-1) on the basis of their theoretical analysis 
and contrail and towing tank experiments concluded that, far downstream 
from the aircraft, buoyancy is the dominant mechanism affecting the growth 
of the wake. They compared their calculations and towing tank data with 
flight data of the visible horizontal width of persistent contrails from 
three different subsonic aircraft. Conti, et al., in reference 8-2, con- 
sidered the possible influences of wind shear, wake collapse, turbulent 
diffusion and buoyancy on dispersion of the wake. For purposes of modeling 
they treated each mechanism separately and by linear superposition they 
concluded that the relative significance of each of the mechanisms was as 
follows t- collapse, 0.33; shear, 0.43; buoyancy, 0.08 and- turbulence, 0.16. 
These calculations were based upon the B— 52 wake width measurements, and 
they postulated that the effect of buoyancy would be relatively insignifi- 
cant when compared to the other dispersion mechanisms. The information 
which is presently available on wake almensions in the final dispersion 
phase is so scant that it becomes a matter of conjecture as to which 
mechanism, if any, is the most significant. 

It is clear from the above discuss chat further experiments are 
required to resolve the problems of buoyant wakes. Flight data of wake 
dimensions together with local stratospheric turbulence and wind conditions 
should help clarify the issue but until these are obtained, we must rely 
on laboratory and wind-tunnel experiments to clarify the problem. 
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9.0 CALCULATIVE EXAMPLE 


9.1 Z n t r oduction — 

in order to show the application of the foregoing theoretical work, 
it will now be applied to an illustrative example of a supersonic aircraft 
flying in the stratosphere.- For illustrative purposes, we will assume an 
airplane of wing span-59 feet flying at an altitude of 70,000 feet at a 
Mach number of 2.75. -It will be assumed that the wing is elliptically 
loaded and the airplane weight is 100,000 pounds. With regard to the 
propulsion system, it is assumed that the engine exits are located on the 
wings with their centerliries_at the middle of the semispans. 

In the ensuing example, we will calculate the wake characteristics 
at the various stages using the theories developed in Sections 4 to 8. 

It is assumed that the atmospheric properties are given by the standard 
atmosphere as listed in reference 9-1. 

9.2 Pure Jet Phase 

A calculative example for supersonic flight has been carried out 
using the theory and correlation curve given in figure 4-2. The following 
initial jet conditions have been assumed: 

M w - 2.75 Mj - 2.53 T Q - 0.308 U w - 0.575 

For future reference the calculated results have been tabulated in 
Table 9-1 -in nondimensional form. This table shows that the length of 
the potential core is about 66.4 nozzle-exit radii. The calculations 
have been carried to a distance farther downstream than the point Where 
entrainment of jets by the vortices has become complete. The column 
headed E^/E^ shows how well the law of the conservation of energy has 
been fulfilled by its closeness to unity. The terms U ave / U oo and 
Tave/^oo Sfive the n> ass average jet velocity and temperature at the various 
axial positions. The ratio V^/U^ shows the magnitude of the radial 
velocity at the edge of the jet. 

For illustrative purposes we show the variations in centerline temp- 
erature and mass average cross-sectional temperature in figure 9-1 sb 
taken from Table 9-1. The centerline temperature is constant until the 
end of the potential core and decreases fairly rapidly thereafter. The 
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average temperature starts to fall immediately because of the high density 
free-stream air being entrained initially. The average temperature is much 
lower than the centerline temperature and decreases fairly slowly after 
the end of the potential core. 

It is possible to determine the local relative concentration of the 
exhaust gases as a fraction of their concentrations at the engine exit 
(i.e., its mass fraction). If we assume that the various gas components 
of -the engine exhaust remain- in the same relative proportions, then the 
relative concentration or mass fraction, c, of the exhaust gases is given 
by equation (4-77) _ The mass fraction, c, has been determined with the 
help of-Table 9-1 to obtain concentration profiles at two axial stations. 
These profiles of relative concentration are shown in figure 9-2 for an 
axial station at the- end of the potential core and an axial station at 
138 exit radii downstream of the exit- 

9.3 Initial Interaction 

The initial interaction equations (5-18) through (5<»21) have been 
integrated downstream to the point of strong interaction between vortex 
and jet using the above initial jet conditions and the following aircraft 
parameters and geometry: 

Jet downward inclination^- 4.9° 

Circulation parameter: K Q ■ 935 ft 2 /sec « 87 meters 2 /sec 

Span: b ■ 59 ft ■ 18 meters 

Nozzle radius^ r* ■ 3 ft ■ 0.915 meter 

n 

with the jet nozzle assumed to be located at half the semispan of the wing. 
The results of these integrations are shown in figures -9-3 and 9-4.. Fig- 
ure 9-3 is a side view of-the vertical positions of the jet and right-wing 
panel vortax centroid, as they proceed downstream, observed from an 
inertial coordinate system attached to the trailing edge of the aircraft. 
The jet closely follows its initial downward slope of 4.9° for approxi- 
mately 12 meters (40 ft) and then, due to both the entrainment of free- 
stream air and the action of crossflow drag, rapidly begins to align itself 
in the free-stream direction. The vortex centroid proceeds downward at 
an initially slow and uniform rate, but begins to move much faster at 
approximately z* - z • r* « 60 meters (200 ft), indicating that the vortex 
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ana jet are interacting strongly at this location. This interaction is 
most clearly observed in figure 9-4 which shows the crossflow plane view 
of the motions. The effect of crossflow drag and ingested free-stream 
air in slowing the initial downward motion of the jet is also evident 
from the relative spacing of the jet trajectory points. While the calcu- 
lation has been carried to z* m 73 meters (240 ft) , it is clear that the 
jet and concentrated vortex have begun to interact strongly somewhat up- 
stream of this location, at approximately z* = 60 meters (200 ft) . 
Consequently, the process of jet break-up and ingestion into the vortex 
recirculation region should be started at or upstream of z* - 60 meters 
for this particular example. 

9.4 Jet Disintegration 

The differential equations (6-1) through (6-4) describing the jet 
disintegration process have been applied to the present example at various 
longitudinal positions at and upstream of the location where the initial 
interaction calculations indicate strong interaction. In these results, 
shown in figures 9-5 (a) through 9-5 (c), the jet has been subdivided into 
five equal-area segments, with initial crossflow plane locations and 
velocities provided by the initial interaction results. In order to faci- 
litate visualization of the ingestion process,, the results are presented 
in vortex coordinates. Figure 9-5 (a) shows the crossflow plane trajec- 
tories of the hot gas segments when the break-up process is initiated at 
z* m 59.6 meters (195 ft). The trajectory points correspond to Az* 
increments of 9.1 meters (30 ft), and the calculations for a particular 
segment are terminated when that segment approaches within 0.3 meter (1 ft) 
of the vortex centerline. These results indicate that all of the hot gas 
segments are captured by the vortex at approximately z* ~ 120 meters 
(400 ft) . 

In order to illustrate the sensitivity of the break-up process to 
the assumed point of its initiation, figures 9-5 (b) and 9-5 (c) have been 
prepared showing the hot gas trajectories when the break-up process is 
initiated at zj « 46 meters (150 ft) and 23 meters (75 ft) , respectively. 
Figure 9— 5(b) shows that when zj — 46 meters, a point clearly upstream of 
the strong interaction location, all of the hot gas is still captured by 
the vortex, with all the hot gas spiralling in toward the core. Fig- 
ure 9-5 (c) displays the analogous results for zj • 23 meters, a point 
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slightly more than one wing span of the trailing edge. These results 
indicate that even if the break-up process is started at this very early 
stage, 60 percent of the gas will be ingested into the vortex recirculation 
region. Consequently, for the parameters and initial conditions of this 
example, the present theory indicates that the vortex pair will capture 
all of the hot jet exhaust within essentially 120 meters (6-7 wing spans) 
downstream of the trailing edge. 

9.5 Equilibrium Vortex Phase 

It is of interest first to look at the characteristics of the 
trailing— Arortices— before ingestion of the hot— jet into its center as a 
guide to evaluating the effects of such ingestion.- First let us determine 
its tangential velocity profile and its temperature distribution. In 
determining its tangential velocity distribution it is necessary to obtain 
an approximate idea of its core diameter and its maximum tangential 
velocity. It has been found in such papers as references 9-2 and 9-3 
that the core diameters are relatively insensitive to angle of attack. 

It has also been found that for a Convair— 990 wing model, the core dia- 
meter is about twice as large as that for a rectangular wing. The loading 
on the wing of the Convair 990 model is probably fairly representative of 
an elliptical span loading. From a measurement of r x - 1.1 inches 
(2.8 cms) for a rectangular wing of 8 feet (2.44 m) span, the scaling by 
a factor of 2 to a 59-foot span wing yields 

r x “ 1.4 feet (0.43 m.) 

With regaro to maximum tangential velocity, Orloff and- Grant (ref. 9-2) 
find that the maximum tangential velocity as- a fraction of V w scales 
linearly with angle of attack for a rectangular wing. This agrees with 
the finding of constant core diameter and a constant nondimensioral 
velocity field in planes perpendicular to the vortex axis with as 

the nondimens ionalizing quantity. Also, the work of Corsiglia, Schwind 
and Chigier (ref. 9-3) show that there is little decay of the peak tangen- 
tial velocity with axial distance up to 31 span lengths behind the wing. 

An estimate of the maximum tangential velocity for both the rectangular 
wing and the Convair 990 model is given by 

- °- e ,9 - 1) 
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This yields a value of v x for the present cases of 

Vl - 271 ft/sec (83 m/aec) 


Thus 

Ki - (1.4) (271) - 379 ft 2 /sec ( 35. 2 m 2 /sec) 

Based on these values and equations (7-1) to (7-3) the nondimens ional 
tangential velocity distribution shown in figure 9-6 has been calculated. 
In the wake region, the extension of the law of the logarithmic region is 
shown since this approximation has been used in the theory of Section 7. 

We now apply equation (7-40) to determine the temperature distribution 
in the trailing vortex before hot gas ingestion. The results are shown 
in figure 9-7. The theory shows that the static temperature reaches a 
minimum at the eye, and is constant in— the—eye. 

Consider now the effect of hot gas ingestion into the vortex. The 
hot gas is completely ingested at 133 exit radii assuming the break-up 
starts at 65 exit radii at the end of the potential core. Between these 
limits ws have the following approximate quantites for the pure jet 


z 

T ave^ T « 


c ave 

n 

(X* - T*) 

v ave »' 


66 

138 

1.53 

1.34- 

1.784 

2.10 

0.205 

0.133 

5.23 ft 
6.3 ft 

208°P (116°C) 
134°F ( 75°C) 

1.6 m. 
1.93 m. 


We now-use these values to estimate the temperatures and concentrations at 
the end of ingestion— During ingestion-the question arises whether the 
rate of mixing with free^stream air- is increased or decreased with respect 
to the mixing rate of the free jet. During ingestion the air tends to 
move into the center and -displace the air in the inner region of the 
vortex. The free-jet calculations show a decrease in average mass frac- 
tion from 0.205 to 0.133 over the length of the ingestion process. The 
effect of the vortex is generally to inhibit entrainment of air into the 
center of the vortex. Accordingly, it la thought that the average mass 
fraction in the vortex center at the end of jet— ingestion will be between 
the two values quoted above. 
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Several other points now need to be made. First the radius of the 
ingested hot gas is now 6.3 feet (1.93 m.) versus a core radius without 
ingestion of about 1.4 feet (0.43 m.) . It is clear that the whole struc- 
ture of the cold vertex can be changed as a result of ingestion. The 
possibility of vortex bursting cannot be dismissed. Additional analysis 
of the vortex with the hot gas core has not been made. However, at 
r ■ 6.3 feet the temperature drop associated with the turbulent vortex 
solution will yield a very small drop in temperature (fig. 9-7) so that 
no substantial cooling of the gas from 134 F above the free-stream tempera- 
ture should occur.- Accordingly, the center of the vortex out to a radius 
of 6.3 feet should have an average temperature approximately 149°F above 
the free stream at the end of ingestion, and all the exhaust gas should be 
trapped in this radius with an average concentration equal to or greater 
than 0.133 of its value at the jet exit. 

One can now ask the question of how fast the hot gas core will cool 
off as the equilibrium vortex phase continues. In order to answer this 
question precisely, we must know how the vortex structure was changed as 
a result of hot gas ingestion. (What changes occur in r 1 and K x ?) 
Unsteady momentum and heat transfer analyses of the vortex are necessary 
to answer this question. Whether unsteady heat transfer effects would be 
significant before the buoyancy phase starts within ten seconds is not 
known. Thereafter the vorticity distribution will be altered by the mass 
entrainment associated with buoyant rise of the recirculating fluid. 

We expect a slow approach of the vortex with a hot core to the thermal 
equilibrium temperature distribution derived in Section 7 through analogy 
with the slow diffusion of vorticity from an equilibrium turbulent vortex. 
The diffusion of heat and of vorticity should have comparable time scales 
for a turbulent vortex. Diffusion of vorticity has been studied in ref- 
erence 9-4, and a numerical example has been carried out. The vorticity 
in the recirculation region can diffuse outward across the enclosing 
streamline or it can diffuse to the vertical plane of symmetry where it 
ia annihilated by symmetry. For the particular example considered in 
reference 9-4, it was found it took more than an hour for the vorticity 
in one half of the recirculation region to decrease by 1/e. It is thus 
probable that the hot gas in the vortex center remains hot and concentrated 
until the onset of the buoyancy phase and possibly for 9ome time thereafter. 
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9.6 Asymptotic Buoyant Phase 


In connection with th^ question of when the asymptotic buoyancy phase 
starts, Pay and Overcamp, ref. 9-9, defin. a characteristic time t c 


t 

c 



(9-2) 


U « air mass density b = airplane span 

r CO 

U* = airplane speed L ® airplane lift 


Their data from contrail measurements, figure 1 of their reference, shows 
that the buoyancy phase, measured by the onset of significant increase in 
contrail width, occurs at a time of about 2t c 

For the supersonic example being considerea herein the characteristic 
time t c is only 4.8 seconds so that the onset of the buoyancy phase 
starts at about 9.6 seconds at a distance about 25,000 feet behind the 
aircraft. (This short elapsed time is another reason that the hot gas 
in the center of the vortex has little time to cool down during the 
equilibrium vortex phase.) 

The initial conditions for the buoyancy calculations include the 
following : 

At = average wake temperature rise due to engine heat 
= initial wake radius (assumed circular) 

p^ ■ initial atmospheric mass density 
i 

T* * atmospheric temperature (constant at 392. 3°R) (218°K) 

a = entrainment coefficient 

p T 

r CO , Co 

p. = - - ■ ■ = initial average mass density of wake 

1 + At 


The value of At can be determined from the following formula baaed 
on the assumptions that all of the engine heat goes into the wake and the 
span ioading is elliptical. 
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(9-3) 


At 


0.288 


r, Wc p 


C D = airplane drag coefficient 
r r •> b 2 /wing reference area 

^ a efficiency of conversion of engine heat into 
propulsion 

c - specific heat of air in mechanical work units per unit 
P mass per degree 

For our example case, the following numbers were used: 

At - 9. 3°F (5. 16°C) 

« 43 feet (13.1 m) 

p - 0.0001389 slugs/ft ? , (h * 70,000 ft) 

"i 

T* - 392. 3°R 
00 

a * 0.3 and 0.9 

Equation (8-24) was integrated using the above input quantities for the 
two entrainment coefficients, and the calculated results are tabulated in 
Table 9-II up to the time the wake reaches a.ts maximum rise.. It is noted 
that the static temperature rise for the present example is much larger 
than the fractional part of a degree usually found for subsonic aircraft. 
This result is a consequence of the greater speed and lower aspect ratio 
as shown by the above formula. Table 9-II shows, however, that the average 
temperature of the wake is almost equal to that of the atmosphere at the 
end of the buoyant rise. 

Inspection of Table 9-II shows that the increase in radius of the 
wake due to entrainment is given closely by ah for the entire range of 
the tables. The rate of entrainment is greater- for the larger entrainment 
coefficient as shown in figure 9-8 (a), but the increase is not proportional 
to a. The vertical rise shown in figure 9-8 (b) is slower for the large 
entrainment coefficient as would be expected because the same initial 
buoyant force must now act on more fluid. The decreased vertical velocity 
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tends to make the rate of entrainment less than proportional to a. The 
results for a - 0.3 are believed to bo representative of a quiescient 
stratosphere, and thus only represent an idealized extreme of the real 
stratosphere. The results for a « 0.9 are considered more representative 
of real stratospheric effects until better data on entrainment coefficients 
for line thermals in the stratosphere are available. 

The question of what the concentration of the exhaust gases is during 
the buoyancy phase cannot be answered at this point. Since the buoyancy 
phase starts after about 10 seconds, it is probable that no significant 
decay of the trailing vorticity^ 1 ) has occurred; and that the hot gases at 
the. center of the vortices have not diffused significantly. With the onset 
of the buoyant rise, a new induced circulation pattern is imposed on the 
recirculating fluid which tends to oppose the trailing vortices. At least 
this is the deduction one would make from the results of Tsar.g (ref. 8-5) 
which show that a horizontal cylinder of denser fluid initially with no 
circulation develops a vortex pair of the opposite sense of the trailing 
vortices as they move under the influence of buoyancy. This generated 
vorticity varies theoretically as t 1 / 3 for constant density external to 
the region (ref. 8-5) . A detailed analysis of the buoyant wake considering 
the initial trailing vortex circulation should be made to study the details 
of the turbulent mixing process. An experimental study of this problem 
would also be extremely valuable. 

It is of interest for chemical studies to know the actual distribution 
of temperatures to which the exhaust gases are subject during the buoyant 
phase. No attempt to determine this distribution was made in the present 
work. Upper and lower bounds on the temperature of the exhaust gases 'an 
however be given. A lower bound on the temperature and mass fraction of 
the exhaust gas can be calculated on the basis of the assumption of homo- 
geneous gases in the recirculation region. An upper b^und would be given 
by assuming that the exhaust gases are frozen in the vortex cores at the 
temperature and mass fraction corresponding to the end of hot gas ingestion 


Vortex bursting would probably not diffuse vorticity from the recircu- 
lation region although the Crow instability could break up the 
recirculating region. 
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Figure 9-1.- Axial distributions of centerline and average temperature. 












Figure 9-5.- Trajectories of hot gas exhaust segments 
during break-up process in vortex coordinates. 










Static temperature drop 

















10.0 CONCLUDING— REMARKS 


Detailed analysis has been made of the ingestion and dispersion of 
engine exhaust products in the. trailing vortex system of a supersonic 
aircraft. While the analysis is not limited to supersonic flight, the 
following concluding remarks are directed principally to this speed regime. 

1. An analysis has been made of the mixing between the hot supersonic 
jets of the ,_supe r s on ic aircraft and the conflowing supersonic stream.- One- 
set of experimental data was found-covering this case, and the eddy- 
viscosity coefficients necessary to carry out supersonic jet calculations 
were determined from an empirical correlation of these data. It was found 
that the rate of mixing between the supersonic jet and the conflowing super- 
sonic stream was an order of magnitude less than would be expected on the 
basis of the subsonic eddy-viscosity, results. As a result, the potential 
core of the supersonic jet is much longer than expected and the exhaust 
gases are subjected to elevated temperatures and concentrations over much 
longer distances than previously expected. 

2. A new theory has been developec for the break-up of the engine 
jets and their ingestion into the center of the trailing vortices from the 
wing. It was found that ingestion started at the end of the potential core 
and that all of the hot gas from the engines was ingested into the trailing 
vortices within two core lengths. On the basis of the free— jet calculations 
the mass fraction of the exhaust gases changed only a small amount over the 
ingestion distance. It is thought that it changed even less during in- 
gestion because of the effect of the vortex in inhibiting entrainment. 

3. It was found that the volume of hot gas ingested into the core of 

the trailing vortices was much greater than the volume of the core of the 

trailing vortices alone. It is probably therefore that the tangential 
velocity distribution of the trailing vortices is substantially changed 
by hot gas ingestion, the vortex thermally insulates the hot gases from 
the stratosphere and prevents further mixing. An analysis of the hot 
trailing vortices was not made nor was the possibility of inducing vortex 
bursting by hot gas entrainment investigated. 

4. It was found that the final buoyancy phase of the trailing 

vortices started only ten seconds after leaving the trailing edge. As 

a consequency, the equilibrium vortex phase was of short duration. 
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5. Buoyancy calculations were carried out based on the assumptions 
that the wake is essentially horizontal with no axial flow. It was found 
that the temperature rise in the wake of the supersonic aircraft is much - 
greater than that for a subsonic transport since temperature rise was 
found to vary directly, as speed squared and inversely as aspect ratio. 

It was found that the calculated results depend on entrainment coefficient 
which is expected- to be dependent for its value on such factors as atmos- 
pheric turbulence, atmospheric winds, and wake turbulence.- Comparison 
between the calculated results for the present supersonic example and the 
correlation of Fay and Overcamp of contrail data for the width of the 
wakes of subsonic aircraft shows good agreement. 

6. In the-buoyant model it was assumed that the recirculation region- 
remains intact during the buoyant entrainment process- Asymptotically the 
solution yields a buoyant mass which rises to a maximum height and closely 
approaches atmospheric temperature. 

7. The buoyancy theory was developed on the basis of constant entrain- 
ment coefficient measured in water tank tests. The value of the actual 
entrainment coefficient for a buoyant wake in the- stratosphere needs to be 
determined by careful flight tests. 

8. One assumption of the theory is that the initial vorticity at the 
beginning of the buoyant phase does not have any significant effect on its 
vertical motion. As a result of entrainment, vorticity is developed which 
tends to destroy the circulation of the trailing vortices. It is probable 
that the theory should -be extended to include these vorticity effects, and 
that -experiments in tanks on line thermals with initial circulation should 
be carried out to see how the entrainment coefficient is affected. 

9. For subsonic speeds it is probably important to develop the 
equilibrium vortex theory with hot gas ingestion since the times involved 
for this phase- are significant. For the supersonic case this matter may 
also be important so that proper initial conditions for the buoyancy phase 
can be determined. 

10. Velocity and temperature profiles have been specified at various 
stages of the wake. If no chemical reaction is assumed to take place, 
concentration distributions of individual checmical species can be deter- 
mined simply from the knowledge of species concentrations at the engine 
exhaust by using the analogy between heat and mass transfer. On the other 
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hand, since the rates of reaction of chemical species are dependent upon 
local temperature and species concentration, the analysis presented in 
this report could.be used to-predict variations of concentration- of species 
under conditions of chemical reaction. In particular, if ambient concen- 
trations of ozone and exit concentrations of NO are measured or assumed, 

local mixture ratios of - NO .id ozone can be determined on the— b asis — that 

ozone is being entrained -into the_exhaust jet* 

Nielsen Engineering & Research, Inc. 

Mountain View, California 
November 1973 
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